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Abstract 

Details  are  provided  of  research  directed  towards  the  development  of  the  computer  pro¬ 
gram  W.AKE  for  predicting  the  internal  waves  produced  by  the  passage  of  a  submarine  in 
a  density  stratified  ocean.  The  research  was  undertaken  within  DSTO  during  the  period 
February-March  1992.  The  report  takes  the  form  of  a  series  of  appendices  (presented  in 
chronological  order),  which  provide  a  record  of  the  progress  of  the  research  during  that 
period.  Results  are  presented  which  indicate  that  a  representative  submarine  produces 
internal  waves  which  have  a  velocity  magnitude  of  about  one  millimeter  per  second. 
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Summary 

In  late  1991.  I  agrt'ed  to  luidenaice  an  investigation  for  DSTO  of  sut.)iuan!!e  wake>  with 
enipliasis  on  the  internal  waves  generated  by  the  passage  of  the  snbnuirme  tliriiugii  a  strat 
ified  ocean.  The  priinary  task  was  to  compute  thtsie  mtenud  waves  fur  a  giveii  sutunarine 
and  this  task  was  completed  in  the  period  Feltruary-March  l')92,  Rt*seaich  a.'.sisrHUce  was 
provided  by  Graham  Furneil  and  Tony  Legg  of  DSTO.  anti  [migTamnung  assistance  bv 
Michael  Carroll  of  Ebor  Pty  Ltd.  Additional  advice  and  guidance  wa-  [irovuitHl  bv  Davai 
CartwTight  and  Dan  Madurasinghe  of  DSTO 

This  report  is  in  the  form  of  a  eolkxtion  of  Appendices,  These  are  presenttsi  m 
clironologrcal  order,  giving  a  record  of  progress  as  the  project  unfolded  Not  ail  of  the 
Appendices  are  of  equal  imixirtance.  and  some  are  pirehtninaiy  or  tentatiw.  in  jian  supiei  - 
sedeti  by  subsequent  work.  All  are  neveret heless  iriciuded.  to  give  an  accurate  impression 
of  the  character  of  the  project. 

■Apfiendix  1  is  the  most  important,  setting  the  style  and  backgrouiiti  for  the  project 
This  Appendix  was  actually  wTitten  before  the  main  jirojetu  commenced  Ajipendix  s 
is  also  very  important,  since  the  project  follows  similar  jiathways  tu  the  Canadnui  effort 
headed  by  Dawson  and  Hughes.  There  are  many  things  riiat  i-ouid  not  fie  done  <ir  done 
sufficiently  thoroughly  in  the  very  short  time  available.  Some  of  thtsse  could  make  good 
topics  for  further  studies.  Appendix  Id  is  a  panial  list  of  rh(>se  topics.  Appendix  Is  i>  a 
very  preliminary  repiut  of  final  results.  It  had  to  be  left  till  the  verv  list,  and  is  presented 
in  the  form  mostly  of  graphs  with  hand-wTitten  annotation.  A  inort'  complete  sununaiy  of 
results  may  be  presented  later. 

The  bottom  line  is  that  a  submarine  at  repiresentative  liepth.-  and  sp(*eds  in  a  rf-pris-eii' 
tative  stratification  makes  internal  waves  of  a  velodtjc  magnitude  of  al'out  oru-  tmUitu<>te! 
per  second.  This  sounds  small,  and  is.  although  (because  of  the  long  wavelength  and  large 
period  -  about  20  minutes  -  of  the  internal  wave*  it  does  correspond  to  actual  [laniclc 
displacements  of  the  order  of  about  one  metre.  These  rtstiits  art'  of  tiie  order  of  magmtude 
of  those  reptrrted  elsewhere.  They  may  or  may  not  be  detectabh'.  at  lea^if  indirt'ctlv  p.g. 
via  bunching  of  capillaries.  Our  program  allows  parametric  studv  of  \'ari<ms  suf>marme 
shapes,  sizes,  depths  and  speeds,  and  of  various  stratificatioiir. 
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Appendix  1;  10  February  1992 


Source  in  Stratified  Fluid 


Some  prelimman’  notes  on  intermil  wave  generation.  a<Tuai{y  prej^ired  hv  EOT  prsoi 
to  commencement  of  project. 


1.  Basic  equations. 

Let  the  fluid  velocity  vector  be  q  =  {1'  —  Oj  Oy.  u  i  where  o  ~  oi.r.  v-  :)  and  if 
u-(T.y.2}  are  small.  That  is.  the  flow  is  a  small  perturbation  to  the  umfurin  stremn  f'  in 
the  direction.  The  above  reprt'sentat ion  implies  that  the  c-i'om[x.iiient  of  the  vorncitv 
is  zero,  but  other  components  are  not:  the  flow  is  not  irrotational.  and  o  is  not  a  velocity 
potential. 

Now  the  continuity  equation  is 


C) X  z  ^yv  ^  ~  b 


1.1! 


and  there  is  a  second  equation  of  motion  for  a  stratified  medium  with  density  p  =  p< :  i.  of 
the  form 


-  Oxxz 


P  . 

—  lOji  --  KtC  r 

P 


(1.2) 


where  ^ 

[I  am  not  too  sure  about  the  reference  statiL'^  of  this  equation  Latn  instrfiorr  In 
fact.  I  now  believe  that  it  i.s  well  documented  and  denved  in  C.-S  Yth'^  iyh5  lx)ok.  .<-> 
bibliography.  Appendix  2.  It  seems  to  be  used  by  people  such  as  Keller,  but  without 
derivation.  I  have  derived  it  myself  in  two  quite  independent  ways,  getting  the  same  answer. 
First  by  linearising  the  Euler  equation  -  straightforward  and  boring.  Set'ond  by  considering 
an  ocean  consisting  of  many  layers  of  uniform  density  fluid  with  standard  Kelvin  linearised 
free  boundary  conditions  at  the  interfaces,  then  letting  the  layer  thicknesses  tend  to  zero. 
Incidentally,  the  so-called  Boussinesq  approximation  seems  to  be  to  drop  the  term  on  the 
right  in  On-  but  I  see  no  reason  to  do  that.': 

Our  task  is  to  solve  these  equations  in  t  <  0  subject  to  the  free  surface  condition 


Oj-j  /CM- =  0  (1.3) 

at  r  0.  and  (for  an  infinitely  deep  ocean)  o.  w  —  0  as  c  -•  -x..  Note  that  the  fr{>e 
surface  condition  f  1.3)  follows  from  ( 1 .2)  by  requiring  the  coefficient  of  the  density  gradient 
to  vanish,  which  is  as  it  must  be.  since  the  free  surface  is  just  a  step  discontinuity  in  density. 

Actually,  we  don't  w^ant  to  solve  (1.1)  evert-where:  instead  we  want  to  allow  a  source 
at  (0. 0.  -h).  Hence  we  should  replcce  the  zero  on  the  right  hand  side  by  a  delta  functicu. 
b(x).b(y).S(2  +  h).  Equivalently,  we  need  .solutions  of  (1.1)  possessing  suitable  singularities 
at  this  point,  see  below. 


2.  Fourier  Decomposition 
Write 


o 


‘^rjy . 


i )  dtid.i 


and 


idad.i 


:2.1 


Then  (with  k‘  =  q*  -  J*  and  using  a  dash  for  duiz).  (Mi  gives  U  '  =  while  1 1.2; 
gives 

ir  ~  -  -£U  )  i2  ;5i 

O' 

Then,  eliminating  we  have  an  ODE  for  IT  =  H  i;,),  iiamelv 


or 


where 


A-'ir  -  H"'  =  L-(ir'  -^u'i 
d 


.pH"d  -  (k‘t)~ap'}\V  =  0 
a  =  Kk~  Q‘  =  K  se<"  d 


i2A) 


(2.0) 


;2.f>l 


with  Q  ■=  k  cos  6/. 

Equation  {2.5)  is  a  straightforward  variable-coefficient  second-order  ODE.  and  is  to 
be  solved  subject  to  suitable  free  surface  and  bottom  conditions  The  former  is  just  the 
Fourier  transform  of  (1.3).  which  leads  to 


\V'  =  a\V  (2.7) 

at  r  =  0.  The  bottom  condition  is  best  expressed  by  assuming  that  the  density  is  uniform 
below  some  level  H.  in  wffiich  case  the  solution  to  (2.5)  for  c  <  -H  is  proportional  to  c*-'. 
Hence  in  that  range 

H"  =  kW  (2.H) 

and  the  appropriate  bottom  condition  is  then  that  (2.8)  also  holds  at  -  =  -  //  for  .solutions 
in  ;  >  -Ef.  Of  course  if  there  were  an  actual  Hat  impeniieable  sea  floor  at  c  =  -  H.  (2.8) 
would  be  replaced  bv  H'  =  0.  but  we  shall  use  (2.8). 

Non-trivial  solution  of  the  ODE  (2.5)  subject  to  both  the  free  surface  condition  (2.7) 
and  the  bottom  condition  (2.8)  is  not  possible  in  general,  but  is  possible  (at  fixed  rr)  for 

special  values  of  the  w'ave  number  k.  i.e.  for  eigenvalues  k  =  h'jia).  j  =  0. 1.  2 .  There 

is  always  at  least  one  eigenvalue  k  —  Koicr)  =  a.  whe'e  the  solution  is  just  proportional  to 

H'  =  e''^  as  is  immediately  apparent  from  the  form  (2.4)  of  the  ODE.  even  if  the  density  is 

an  arbitrary  non-constant  function  of  z.  If  p  is  constant,  this  is  the  only  solution,  and  i.s  the 
usual  Kelvin  surface  wave.  If  p  is  not  constant,  there  may  be  more  solutions  j  =  1,2 . 


iuid  these  are  what  are  called  internal  waves  Tlie  retatiunslup  A'  =  K,{(7)  is  a  form  of  tlie 
dispersion  relationship  for  the  j'th  internal  wave 


3.  Solution  Method  for  Sources 

Although  there  are  main-  direct  methods  for  solving  the  above  eigenralue  proiuem  if 
all  that  we  needed  was  the  dispersion  relation  A-  =  A'.icr).  the  following  indirect  inetli(>d  i> 
preferable  if  we  need  solutions  for  soiirces. 

Suppose  il'  =  and  li'  =  iro(;i  are  two  separate  solutions  of  i2  d<‘hned  a.' 

follows.  For  U'l-  we  satish-  the  bottom  boundary  condition  ''i.M  but  not  'he  free  surface 
boundary  condition  (2,7!.  and  normalise  (arbitrmilyi  the  value  of  If'i  at  tjie  bottuni 
obtaining  an  initial  value  problem  to  be  solveii  upward  in  c  >  ~H.  stmting  with  the 
initial  conditions 

U-;,-//;  =A>-^^' 

at  ;  =  -H. 

Similarly,  lib  satisfies  the  free  surface  condition  (2.7'  but  not  the  bottom  coiKlition 
{2.8).  and  is  normalised  (arbitrarily)  to  the  starting  value  -1  at  :  =  0,  Then  we  soivt' 
downward  for  H  2{r )  as  an  initial  value  problem  in  <  0.  starting  with  the  initial  conditions 

UbfO)  =  -  1 

id.  2' 

ir,  (0)  =  -  tr 


at  ;  =  0. 

Each  of  U'l  and  Ub  always  exists  for  any  choice  of  h.rr.  for  reasonalde  p(:j,  and  is 
readily  computed  bj-  any  standard  ODE-solver  numerical  package.  Unless  k  and  rr  are 
connected  by  the  dispersion  relation.  If)  and  lU^  are  linearly  independent  of  each  other, 
but  if  it  happens  that  k  =  Kjiu).  they  must  be  proportional  to  each  other,  since  then  and 
only  then  can  both  boundary  conditions  be  satisfied  by  the  same  function 

The  Wronskian  IU1U2  -  llbU','  must  be  inversely  proportional  to  the  density  pi:),  a 
standard  property  of  variable-coefficient  equations  like  (2  -S).  so  we  can  wTitc  for  all  c 

H-i(;)li^(;)  -  »2(;)IU;(:)  =  Dp{-h)^p(:)  (3.;U 


for  .some  constant  D  =  Dik.n).  the  value  of  the  Wroaskian  at  r  =  -h.  Now  consider  the 
following  discontinuous  solution  of  (2.o).  namely: 


»o(rj 


-ir,'(-/>)Ub(c)  if  ~h  <  i  <  0. 


(3.4) 


[Note  that  lid.  like  li'i  and  U^.  depends  implicitly  on  the  parameters  k  and  rr  as  well  a'- 
the  coordinate  2.  and  will  be  displayed  as  lld(r:  k-rr)  whenever  it  is  necessarv  to  indicate 
that  dependence.]  The  function  iid(-)  satisfies  both  boundary  conditions  (2.7)  and  (2.8). 
by  construction.  Its  derivative  is  continuous  across  2  =  -h.  whereas  its  \-alue  jumps  by 
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-1-  0)  -  \\\)(-h  -  0)  =  D.  Houce  M'(;)  =  Uo(c)  D  is  the  required  diseouTuimjus 
solution  for  a  unit  source,  which  must  have  a  unit -magnitude  jump  across  ;  =  -h 

The  important  point  is  that  D  =  0  when  k  ~  K,[iy).  Tliat  is.  linem  mdejiendeace 
fails  when  the  Wronskian  vanishes.  Hence  D{k.o)  ~  0  Is  an  implicit  form  of  the  disjier.Miui 
relation,  and  that  relation  can  he  detenmnetl  by  numencal  solution  for  U'l  and  U';,  and 
hence  for  D\k.o).  The  solution  U'o  is  •a'ell  defined  and  bounded  for  all  k.cr  valiu's.  Hence 
the  source  solution  U  o  D  is  singular,  with  a  {lole  wherever  D  =  1).  i.e.  wherever  k  K.  icr ) 

'Figuiv  1.1  i.s  a  sampk  plot  of  the  Wron.skian  D(k.rr)  istx  ixL.H  <it  :  ■---  h  r.oj.'.  k 

at  various  a.  Note  the  value  where  D  —  0.  See  Fipuir  6.1  for  the  vont  sporidnit}  dt.spt  r.‘^ton 
relation  curres  k  =  Kj{a).i 


4.  Pole  Avoidance  and  Free  Waves 

The  .solution  now  found  for  the  vertical  velocity  w  is  given  by  the  Fourier  integTai 
i  l.a'(.  which  after  a  change  to  polar  wave  numbers  by  o  =  A*cos«.  .i  —  k<m6  tiecomes 


w\.r.  y.  z) 


_  I  I 


Uni  r:  k.  s  >ec'  6 ; 
Di  k.  X  sec*  t>) 


kdkdh 


However,  since  D  =  0  at  A-  =  Kj{x>ev~  6).  we  niu,-=t  di.ston  the  path  of  A'-integrarion  to 

avoid  the  poles  where  D  =  0.  This  proccd'ire  i?  the  .'ame  f(>r  all  jtoles  j  ~  0. 1.2.  .'5 . and 

is  well  established  for  the  Kelvin  surface  wave  pattern  j  =  0.  Nimieiv,  if  cos  6  >  0.  <list(,)n 
the  path  above  the  pole,  and  if  cosB  <  0.  distort  it  below  the  jiole. 

The  reason  for  this  choice  is  clear  if  we  consider  rotation  of  the  path  of  A'-integration 
through  r90'.  If  x  <  0.  and  this  rotation  is  done  tlirough  -90"  when  cos «;'().  the  {lole 
is  not  cros.sed  during  the  process  of  this  rotation,  and  the  resulting  integral  along  tlie 
imaginary  Ar-axis  contains  an  exponentially  decaying  factor.  Hence  w  tends  to  2t>ro  rapidly 
as  X  —  -  X.  and  there  are  no  waves  upstream,  as  required  by  the  radiation  condition. 

Once  this  decision  about  pole  avoidance  is  made,  omr  task  of  determining  the  flow 
from  a  source  is  in  principle  completed,  and  it  •‘only  "  remains  to  evaluate  the  double 
integral  (4.1).  with  U'o  and  D  known  by  solution  of  the  ODE  (2.5)  (for  all  rallies  of  k  and 
a  =  Ksec^ft).  This  is  a  raanmioth  task,  bearing  in  mind  that  only  recently  ha-  it  even 
been  considered  feasible  to  do  routine  computations  for  unstratihed  fluids,  where  the  ODE 
part  of  the  task  is  eliminated. 

However,  the  main  far-field  (x  —  -f-oc)  contribution  is  from  the  pole?.  This  contribu¬ 
tion  can  be  estimated  by  rotating  the  path  of  A'-integration  in  the  opposite  direction  from 
that  for  X  <  0.  In  that  case,  one  passes  across  the  poles,  so  picking  up  a  contribution 
from  their  residues,  before  arriving  at  an  integral  on  the  imaginary  A'-axis  which  again 
tends  to  zero  rapidly,  as  x  —  -i-oc.  Hence  the  dominant  terms  m  w  as  x  -*-x  iire  the 
contributions  from  the  poles,  namely 


U'f  (x.  p.z  ) 


^  j  2771  sgiifcos 0 j  ^  e'*'"  ^ 


Dk  {k.x  .sec^  B) 


where  k  =  Kj{ksoc'~B). 
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Xore  that  the  denoniinator  in  (4.2)  i>  the  partial  /.’-derivative  of  the  Wnnb-kian 
Dik.ij).  evaluated  at  the  point  k  =  h'A.f^)  where  D  -  0.  and  thi>  (luantity  nm<t  In- 
evaluated  nuniencally  with  prevision  as  part  of  the  over-all  voiuputatioiial  t;isk.  Tlie  above 
expression  van  be  simplified  slightly  to  a  real  integral  on  lO.  r  2i.  iianieiy 

ii'f' =  ——  >  /  sin(A';.rvos//)vosi  Ayysinf/i-^  A’,  (/ft  4.4^ 

~  ^  Jo 

The  expressioiL'i  i4.2)  or  (4.'V>  are  single  integrals  wirii  resjievt  to  e.  and  obvuiu.'iy 
an  order  of  magnitude  easier  to  evaluate  than  the  diiiihle  inti'gTal  s4.1  i  In  the  slup-waw 
example.  fVf  is  usually  vailed  the  ■'free-wave"  vontribution.  and  we  retain  that  rerniinolngy 
here.  Although  free  waves  are  only  a  far-held  approximation,  the  ship-wa.e  exjienenve  v 
that  they  approximate  the  full  flow  held  verv  well,  as  viost'  as  one  or  two  shifi  lengtiis 
astern. 


5.  Stationary  Phase 

Although  I  advovate  evaluation  of  the  integral  (4.4)  without  further  appro.xjmaTion. 
it  is  possible  to  make  the  usual  stationary  phase  approximation  for  larg(’  r  -  '.r*  —  ,v' ' 

observing  that  (4.3)  is  the  sum  of  real  pans  of  integrals  of  the  form 

j  dH  cTl 

where  the  amplitude  function  F(H)  is  well  behaved,  and  the  phase  function  is 

«<(//)  =  A'(tf)eos(f/  -  *.),  ..',2,. 

where  --  =  arctan(,!//.r)  is  the  polar  angle  in  the  horizontal  plane,  and  K{d)  ==  K,  (s  <e<-'-’  ft  i 
For  large  r.  the  main  contribution  to  the  integral  (5.1;  is  from  the  neighbouriiood  of 
stationary  values  of  6.  namely  those  where  't'fff)  =  0,  whicli  satisfy 

tanf^  - -■)  =  A"(^)(  Ai{i^)  '*.3) 

For  most  reasonable  K(0)  there  are  two  roots  B  of  (5.3)  for  each  Ipsj-  than  a  verram  upper 
bound  *(().  and  none  above  that  value.  The  twf)  roots  correspond  to  transverse  (smalhu  ft) 
and  diverging  (larger  0)  waves  which  are  observed  for  ;-c  <  -  and  there  are  no  waves  for 
!'i  >  '•'o.  The  waves  tends  to  be  greatest  in  magnitude  near  >■  =  Deternunatiori  of 
fhe  value  of  is  one  of  the  important  tasks,  and  this  value  is  different  for  each  internal 
wave  mode  j.  If  we  re-write  (5.3)  as 

-■  =  H  -  axvX’axiiK' (B) !  K(B))  (5.4) 

we  see  that  the  upper  bound  for  h  '  luust  occur  when  d^'  jdB  —  0.  which  leads  to  a  condition 
involving  the  second  derivative  K"{0)  that  is  not  worth  wTiririg  down.  For  general  K(B) 
it  is  best  to  simply  compute  -  .{d)  from  (5.4)  and  note  its  maxima  or  mimrna. 


For  example,  the  Kelvin  surface  wave  j  --  0  ha.-;  AV  ~ 
becomes 


'  —  H  -  arctamCtanfi) 


<T.  si.i  KW)  ■■  X  sec- W.  aiul  '  a  4? 


which  is  negative  for  positive  with  a  mmnutim  given  by  ---  -  =  mcraml  '1\  2  ^ 

19.5'  when  B  =  arcTan(I  v  2)  36'.  This  value  of  ".i  is  the  famou.'  Kelvin  ^hlp  wave- 

angle.  and  apphes  at  any  spewed  for  any  moving  object  at  am-  suiimergeiice,  prudui'ing 
waves  on  the  surface  of  an  infinitely  deep  Huid.  even  if  the  fluid  lietisity  i.-  non-consfanr , 
But  the  internal  waves  j  =  1.  2.  3. . . .  will  have  a  different  and  there  ,<t>eni>  to  be 
evidence  that  it  is  smaller,  i.e.  that  the  internal  wave  wakes  iu'e  mirrower  tiiaii  Kehin 
wakes. 


6.  Constant  Density 

The  special  case  of  constant  density  is  worth  giving  in  full  Then  there  is  oniv  the  oik' 
mode  j  =  0.  and  we  can  make  the  followine  idem ihcat ion.' 

=  ‘6.1) 


s>'.r}  =  -■  coshiA'c  i  -  -  sinli(  A';' 

K 

D{  Ac  (7 )  =  A'  -  (7 


llo! r;  A',  cri 


A'  -  rr 


,.A*  Z—h 


k  -  O' 


.•<gnl  c  -  lot 


and 


sin{«Tsec6'lcos(/r,i/.sec'^i'sin6')e'' ■  ^ 


which  is  well  known. 

FiguTP  1.2  sh<)u:s  computations  from  f6.5j  by  D  Madumsmohe., 


>12: 
i6.3i 
6.4  t 

'6,5  i 
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Appendix  3;  18  February  1992 

Dispersion  Relations 


I'Note:  In  part  improved  by  later  work,  .sw  .Apjwndtx  (>  to  follow 

Our  concern  here  is  very  much  with  dispersion  relations  for  internal  waves,  nanielv 
relations  between  wave  speed  c  and  wave  nuiniier  k.  For  i-onvenience,  we  use  instead  <>;  < 
a  quantity  proportional  to  its  reciprocal  square,  naniely 

(T  —  g,  r* 

where  g  is  gravity.  Then  we  need  a  connection  between  k  and  n  e.g.  k  --  Ko-^)  Oiie‘ 
of  our  first  tasks  is  to  compute  this  relation  for  a  given  density  distribution  R'cau.-.e 
this  numerical  task  obscures  some  qualitative  features,  this  report  sunmiariscs  some  ^uch 
features. 

The  empirical  dispersion  relation  attributed  to  Phillips  11977)  by  Tulin  and  Miluh 
(1990)  can  be  wTitten 

Ao 

— -n  =  k \coxh(k  H —  1  —  /re 
P  ■ 

This  contains  as  its  leading  term  the  finite  depth  dispersion  relation 

(7  —  kcoxhlkH) 


for  water  waves  in  a  fluid  of  constant  depth  H.  but  modified  ^  ;th  gravity  redui'ed  in  the 
ratio  of  the  relative  density  change  across  the  thermocline.  The  parameters  are  the  depth 
H  from  the  free  surface  to  the  beginning  of  the  thermoi  line  and  the  thickiu^s  (  of  the 
thermocline. 

[Figure  3.1  is  a  .‘Sample  of  the  k  -  a  relationship.  It  is  not  ba.<ted  on  Phillips  formula, 
but  IS  similar.. 

The  small  and  large-ir  limits  are  of  interest.  For  small  k.  wt  have 


a  — 


k^oik‘y 


so  if  <7i  is  the  value  of  a  for  the  first  internal  wave  at  /r  =  0  then  we  intimate 


-1 


and  then  for  o  >  os  (but  close  to  it)  we  have 


k  —  —fa  -  17] )  -  0{(T  -  Oi  i* 

P 

This  apphes  only  to  the  first  internal  wave  mode.  Phillips  doesn't  speciT-  it.  but  I  believe 
for  the  second  mode,  =  9(Ti  approximately,  and  more  generally  that  the  th  mode 
starts  at  t7  =  (Tj  where  approximately  =  {2)  -t-  l)^<7i.  That  is.  I  expect  that 

k  —  Cj{u  -  a j)  +  0(a  -  Oj  }* 
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for  some  Cj . 

The  large-A-  limit  is 

!>( 

That  is.  (7  and  k  grow  large  together,  with  k  \‘arying  tu;  the  squtirt*  root  of  .7.  Tla*  coefhi  lont 
of  proportionality  \~aries  like  the  squiire  riKtt  of  the  density  difference,  and  also  like  the 
inverse  square  root  of  the  thermorhne  thickness  t.  but  does  nor  depend  on  the  tiiermo<-lm<' 
depth  H. 

I  am  rather  interested  in  corrections  to  this,  espwiaily  those  which  (hx-nininate  be¬ 
tween  modes.  I  am  inclined  to  think  that  rr’'  -  needs  to  be  replaced  by  le  -  7*  -  for 

some  cT*  which  is  not  the  same  as  aj  though  it  must  incretise  (like  <7j  dot's  a-^  tiit'  stjuiire 
of  the  mode  number  j.  Also.  I  think  that  in  general  there  should  be  an  additive  term,  .vi 
the  general  large-A  expansion  should  be 

k  —  Ci/T  -  ■  -  kj 

where  C  is  a  universal  coefficient  independent  of  mode  number,  and  a*  and  A,  tiefiend  011 
mode  number. 

One  cannot  rely  on  the  Philhps  formula  for  this  type  of  additional  information  since 
it  assumes  a  special  form  for  the  density  distribution  a<  evidenceti  from  the  fact  that 
the  results  depend  on  parameters  like  H  and  e  which  are  only  really  meaningful  in  the 
context  of  that  distribution.  So  a  useful  piece  of  anahuic  or  .senn-numerical  work  would 
be  filling  in  some  of  these  details  for  a  general  density  profile.  Thar  is.  what  parameters 
of  the  profile  do  quantities  like  <r^.  a‘.  kj.  Cj.  C.  etc.  depend  upon'.’ 

The  Philhps  empirical  formula  is  an  example  of  what  one  can  do  with  a  so-<‘alled 
Boussinesq  approximation.  This  makes  use  of  the  fact  tiiat  the  densitc-  variati()n>  are 
small  in  absolute  terms.  It  only  takes  them  into  account  when  multiplied  by  the  gTavir\- 
parameter  g.  Then  one  can  also  use  a  large-cr  approximation.  In  effect,  one  lets  rr  tend  to 
infinity  and  u  tend  to  zero,  while  retaining  the  product  a.i/.  There  are  some  simplification^ 
to  the  general  ODE  problem,  though  not  too  many.  I  think  there  is  a  lot  could  be  done  of 
this  nature. 

Note  that  the  general  Philhps  formula  has  this  Boussinesq  charai  ter  My  disfuission 
above  pretended  that  was  arbitrary,  but  in  tact  the  Philhps  formula  is  only  valid  if 
it  is  small,  in  which  case  a  is  always  large,  of  the  order  of  the  inverse  square  rtxit  of  Ap. 
whereas  A  is  not  necessarily  small,  ranging  from  0  to  in  general. 
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Appendix  4:  27  February  1992 

Program  A 

The  program  for  submarine  iurerual  wave  wakes  is  at  a  poiiir  of  development  w.hcre 
we  need  to  standardise  it  in  some  way.  For  detimteiu'ss.  let  ii'  l  all  tlie  present  vt'rston 
Program  A. 

Program  A  takes  as  input  a  file  containing  valtit>s  of  the  density  pi ; }  It  ue«‘d.'.  no  otiier 
input.  It  solves  the  ODE  (2.5)  of  Appendix  1,  using  the  Ruiige-Kutta-.N'ysrrom  rnetiiod 
•Vofe.  Ail  equation  numbers  relate  to  that  Appendix. 

The  solution  is  done  only  once  for  each  k.  o.  namely  yielding  the  U'l  i  r  ;  innctuiii  winch 
satisfies  initial  conditions  (3.1)  at  the  bottom  r  =  -H  When  we  reach  the  top  :  =  (i,  wt' 
immediately  evaluate 

Dik.a'i  =  -  (TU’ifO’. 

which  happens  to  be  the  same  as  that  defined  by  ecpiation  (.3.3!  when  r  =  d  and  h  ~  (i 

Thus  we  have  a  function  subroutine  for  Dik.n).  This  is  now  used  oiilv  to  find  its 
zeros.  That  is.  we  use  the  secant  method  to  solve  the  equation  Dik.at  ~  0  Tins  i  iui 
be  done  either  by  finding  the  value  of  a  at  fixed  k  or  (as  m  the  report  )  fry  finding  the 
value  oi  k  ~  Kia)  at  fixed  a.  Each  internal  wave  is  represented  by  a  separate  function 
J  =  0. 1. 2. . . .  and  we  can  find  each  .solution  separately  numericailv.  However,  wi- 
probably  are  only  interested  in  j  ~  1. 

Once  K{a)  is  determined,  the  only  thing  the  program  pres<uitly  does  is  to  plot  n. 

Program  structure 

Presently  a  collection  of  modules.  Xamelv  a  main  program  SEARCH'  and  subpro¬ 
grams  "DET-.-FIND  INTERVAL"  ■RTSEC", "DENSITY". 

These  are  linked  by  files  called  MAKEFILE.  COMP  etc.  which  do  some  things  Mike 
should  write  up  (2-3  paras).  There  is  also  a  library  called  "K  S  GRAPH”  and  some  other 
bits  and  pieces. 

The  irafKirtant  subprogram  is  “DET".  This  does  the  ODE  solving,  by  calling  a  sul>- 
routine  called  "NYST".  which  has  a  subroutine  called  RHSF"  for  the  right  hand  side  of 
the  2nd  order  ODE. 

Program  B-f 

Over  the  next  6  days  w-e  need  to  extend  this  program.  First  (and  this  is  done  already) 
we  need  to  add  determination  of  li  2(3).  This  is  done  by  duplicating  and  extending  "DET". 
Keep  "DET"  almost  as  it  is.  it  is  still  needed.  Write  a  new  subroutine  "WTUNCTS"  which 
computes  both  Hi(2)  and  Il2(~)-  This  is  all  it  dt^es.  and  it  outputs  these  functions  and 
their  derivatives.  It  could  also  compute  the  Wronskian  UMniUhTr)  -  UofrlU','':)  if  we 
wanted  to  do  that  in  this  subroutine,  though  perhaps  that  is  better  left  for  later 

Importantly,  all  the  above  is  done  with  only  the  density  p  as  input,  togfUher  with  the 
parameters  k.a.  (And  some  numerical  or  tolerance  parameters'.').  We  do  not  yet  have  to 
specify  anything  about  the  submarine,  in  particular  not  about  its  depth  h  of  submersion. 


Now  we  specih'  h.  (It  ha.s  to  be  one  of  the  ;  values  for  which  we  iiave  done  the 
computations  of  11'!, ir-.)}.  Then  we  use  (3.4)  to  set  up  Uoi.:).  That  should  be  a  subroutine 
all  of  its  own. 

Now  turn  to  (4.3).  Thar  contains  a  function 

A{z.»)  =  U'oulA'  Da 

which  we  must  compute.  The  0  dependence  conies  about  be<.‘ause  ail  3  cpiantities  tin  the 
right  depend  on  k  and,' or  ct  =  h  sec- 

Here  Hd  is  just  computed,  and  Kirr)  is  found  by  Program  A.  and  Da  is  the  paitial 
derivative  of  D  with  respect  to  k  at  fixed  a.  evaluated  at  k  —  A'(/f  i.  I  think  the  subprogrcun 
"DET"  can  be  modified  to  yield  Da.  since  differences  of  D  are  already  used  in  the  st'cant 
method.  Note  that  we  only  need  Ud  and  Da  with  k  ~  Kia).  not  with  a  genera]  k  and  i7. 

Simultaneously  (Graham)  we  need  to  push  the  program  forward  on  a  couple  of  other 
fronts,  namely  to  go  from  the  K{a)  determination  (subprogram  'DET"  S  toward  kinematic 
aspects  of  the  wave  pattern  (i.e.  those  not  needing  to  know  the  actual  wave-generating 
efficiency  of  the  submarine)  One  immediate  objective  is  the  Kelvin  stationary  phase  angle, 
thus  the  versu  6  information  of  equation  (-5.4j.  Other  interesting  matters  are  generali¬ 
sations  of  what  is  in  the  Tuck,  Collins  and  Wells  (1971)  Journal  of  Ship  Rtsst^arch  paper. 
In  particular,  generalisations  of  the  frequencies  (Fig  6)  obser\'ed  by  Fourier  transforming  a 
one-dimensional  cut,  and  the  “ridge''  or  expected  location  of  peaks  m  the  double  Fourier 
transform  of  a  two-dimensional  record  (Fig  7). 
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Appendix  5:  2  March  1992. 

EOT  Summary 


Here  is  my  suimnary  of  the  submaiiue  wake  program  so  fm. 

My  fall-time  involvement  began  10  February  1992  aiul  is  due  to  end  27  Man  h  1992. 
The  initial  team  with  full-time  involvement  consisted  of  Tuck.  Furnell.  Legg  and  Maduras 
inghe.  from  DSTO.  plus  Carroll  on  contract  from  Ebor  as  programmer.  Part  time  imoht'- 
ment  or  overseeing  role  is  provided  by  Cartwright.  Haack.  Webster.  .Marwood.  etc 

Division  of  responsibiliry  initially  seenietl  to  be  Tuck  for  general  ideas.  Furnell  for 
long-term  DSTO  cariying  on  of  tht'se  ideas  after  27  March  Sand  hence  close  hason  witli 
Tuck  before  27  M-urh).  Legg  for  input  and  (espe<-ially  grapliici  output  considerarions 
Madurasinghe  for  the  numerical  integration  end  phase  of  the  numerical  work.  Carroll  duo 
programming  directly  as  required  bv  Tuck. 

The  team  immediately  seemed  too  large  to  me.  Hence  at  my  suggt'st  ion.  the  immediate 
working  team  was  reduced  to  Tuck.  Furnell.  Carroll,  at  least  fur  the  first  few  wwks  This 
report  summarises  progress  by  that  team. 

For  the  first  two  weeks,  the  team  was  mostly  feeling  its  wa\’  and  getting  to  know 
each  other.  I  arranged  to  transfer  some  veiy  crude  computer  programs  here,  and  Carroll 
started  to  modify  these.  He  was  at  first  looking  at  ways  to  improve  the  ODE  solving 
aspect,  using  non-umform  grids,  but  this  proved  a  non-useful  approach,  espxviaily  since 
we  intend  to  use  a  uniform  density  grid,  so  was  abandoimed  in  favour  of  retaining  the 
Runge-Kutta-Nystrom  method  in  my  original  program. 

I  also  brought  some  Macintosh  programs  which  used  graphic  input  and  output  tech¬ 
niques.  These  were  of  interest  mainly  to  indicate  to  all  members  of  the  team  what  sort  of 
thing  we  wanted 

.A.t  this  stage.  Furnell  was  mostly  concerned  with  parallel  development  of  anal>-"ic 
models  for  step-wise  exponential  density,  ma;  y  to  use  as  che(.’ks  on  the  general  program. 

On  19  February.  I  gave  a  seminar.  After  this,  the  focus  of  the  work  was  clearer,  and 
the  reduced  team  set  out  with  more  specific  tasks.  Cmroll  cleaned  up  the  program 
it  presently  stood,  and  this  was  designated  "Program  A"  on  27  February.  This  program 
takes  a  given  density  profile  as  input  data,  and  computes  the  dispersion  relation  k  =  Kin) 
of  the  internal  waves  for  that  stratification.  It  does  so  by  solving  the  ODE  (2.5)  for  ICif;) 
(numbers  refer  to  the  Tuck  prehminaty  note  "Source  in  Stratified  Fluid".  Appendix  1). 
Note  that  it  does  not  .solve  for  Ub  (see  (3.2))  yet.  The  Wroaskian  D{k.n)  i.s  computable 
at  .:  =  0  (for  submarine  depth  h  =  0)  (see  (3.3))  without  the  need  to  compute  H't.  Then 
the  equation  D{k.n)  =  0  is  solved  for  k  at  each  fixed  a  by  the  secant  rule. 

This  dispersion  relation  is  of  independent  interest.  WV  have  devoted  some  time  to  un¬ 
derstanding  the  nature  of  the  k  versas  n  cur\-ps  for  ideali.sed  and  actual  strarificatioas.  For 
realistic  actual  stratifications,  the  n  values  are  very  large,  typically  measured  in  thou.sand.s. 
when  k  is  of  the  order  of  unity.  The  first  internal  wave  starts  with  A-  =  0  at  some  n  =  n-[ 
value  of  this  large  size,  then  k  increases  monotomcally  wuth  n.  eventually  becoming  large 
when  n  is  large,  asymptotically  like  the  .square  root  of  n.  Similarly,  the  second  internal 
wave  starts  with  an  even  higher  value  of  rr  =  <72  9ni.  and  has  a  k  value  always  smaller 
than  that  of  the  first  wave.  etc. 


There  are  many  consequenees  for  the  final  form  of  the  internal  wavts  of  this  character 
to  the  dispersion  relation,  and  Fiirnell  is  pursuing  some  of  tiiese  for  realistic  ami  idealist^d 
stratifications.  These  include  the  nature  of  the  Kehin  wake  pattern,  in  particular  its 
angle  (narrow  widces.  mostly  diverging  wave's,  etc.),  and  sonit'  aspects  of  Fourier-analvsed 
detection  as  in  the  Tuck.  Collins.  Wells  (1971)  pajK'r.  Note  tliat  this  work  can  jiroceed 
even  without  a  knowledge  of  the  wave-amplitude  generation  capai>ilit\  of  tin  particular 
body  (submarine)  making  the  waves. 

The  next  phase  of  development  of  the  program  to  determme  actual  wavis  arn[)htudes 
is  to  de\-elop  Program  A  further,  to  compute  H'j.  11',,  (st'e  i3.4ii  and  a  (luanrirv  which 
appears  in  (4.3)  but  is  not  given  a  specific  symbol  there,  namely 

.4ir:(T)  =  A-irntr:A'-fT)  DeiA'.o’t 

Actually  as  written,  this  quantity  apftears  to  dt'pend  on  k  too.  but  we  only  wruit  its  value 
when  k  satisfies  the  dispersion  relation  k  =  Kjlrr],  This  .4  is  the  amphrude  of  the  wave 
generated  by  a  unit  so'ure  ar  the  current  value  of  the  wave-siw't'd  jiarameter  .t.  Note  tiiat 
if  the  actual  wave-speed  is  c  them  n  =  y  e'  Also  note  that  .-ince  c  -•  where  f 

is  the  submarine  speed  at  an  angle  B  to  the  wave  direction,  this  confirms  '2. (in  iianielv 
(T  =  Ksec-^.  with  N  ==  Ultimately  we  shall  consider  .4  A\B)  and  integrate  ^4  .'U 

with  respect  to  B.  so  summing  up  intenial  wave's  of  all  dirt'crions 

Note  that  there  is  a  separate  contribution  to  the  integral  (4, .3)  for  each  internal  wave 
mode  j.  Program  A  simply  computes  as  many  internal  waves  j  a'  w('  specih'.  Bur  we  .'hail 
mostly  only  be  concerned  with  one  or  at  most  two. 

An  important  ingredient  in  the  determination  of  .4  is  the  partial  derivative  Dk  of  the 
Wfonskian  with  respect  to  k  at  fixed  <7.  This  is  the  slope  of  the  D  versa-  k  plor  a-  it 
crosses  D  -  0.  All  three  ingredients  k.  U  o.  Dk  of  .4  are  pre'sently  being  computwi  and 
computation  of  .4  is  being  done,  though  at  time  of  waiting  insufficient  ch«x-ks  haxa-  btx'ii 
made  on  it. 

Next  step  (assuming  that  A  for  a  source  is  computable)  is  to  move'  from  a  .source 
to  a  submarine.  First  we  do  a  Rankinc  ovoid,  namely  a  somree-sink  pair  separated  by  a 
distance  of  the  order  of  the  submarine  length.  When  that  length  L  is  large  compared  t(> 
the  maximum  radius  R.  the  required  source  strength  i.s  (according  to  slender  body  theory  i 
UzR^.  That  is.  we  take  the  unit  source  result,  already  coinputf'd.  multiply  it  by  i'  times 
the  maximum  cross  section  area  r/?‘.  and  subtract  the  same  thing  shifted  bv  a  distance 
i  aft. 

In  fact,  doing  this  explicitly  in  (4.3)  gives  a  result  where  the  “sin''  is  replaepii  by  a 
■‘cos',  and  the  "A"  for  the  unit  soimce  is  multiplit'ii  by 

2sin(^K(i7)Lrosf^) 

the  first  square  bracket  being  the  source  strength,  and  the  second  the  source-sink  separation 
factor.  An  interesting  formal  special  case  is  the  limit  as  L  0  and  /?  —  x  with  R'^L 
bounded,  which  yields  a  submerged  sphere,  with  the  'sin''  of  the  second  square  bracket 
replaced  by  its  argument.  This  may  be  a  useful  test  case  to  compute  first. 


Appendix  6:  5  March  1992 

Short  wave  limit  of  dispersion  relation 

Make  the  Boussinesq  approximation  atid  wTite  u  =  -p'  p.  Tlum  the  UDE  v2.->)  uf 
Appendix  1  becomes 

li'"-  ir  =  0 

There  is  no  hope  of  a  solution  satish’ing  rea'ionaV>le  top  and  botttmi  eonduions  \mless  the 
coefficient  in  square  brackets  above  takes  at  least  some  positive  values.  After  all.  if  it 
wholly  negative  the  solutions  are  exponential- like,  and  can  onh-  have  one  zero,  whereas 
when  the  coefficient  is  positive,  the  solutions  are  sinu.soidal-like,  and  can  have  lots  of  zeros. 
T  beheve  that  there  is  a  theorem  to  that  effect  for  ODEs.  imt  it  is  mtuitivelv  obvious 
anyway,; 

Xow  if  we  let  A-  —  oc  at  fixed  a.  eventuallv  the  coefficient  must  l.>ecome  wholh'  negative 
for  ail  Hence  there  will  be  no  solution  satisfi-iug  the  boundary  <-onduions.  'i  herefore. 
fT  must  increase  with  k  at  such  a  rate  as  to  keep  at  least  some  positive  values  for  the 
coefficient.  This  will  happen  last  (as  we  increast"  k  at  fixe<l  rr)  at  those  :  vaJui's  whert'  lO : ) 
is  greatest . 

Hence  suppose  the  maximum  value  of  n(:)  is  i/j^  and  that  it  occurs  at  r  =  Then 
it  is  clear  that  as  k  and  <7  become  large  together. 

ttx/r,,  -  k‘  —  0 

or 

k  —  v'^^m 

Thus  the  large-A:  limit  of  the  dispersion  relation  depends  only  on  the  behaviour  of  the 
stratification  at  the  depth  where  the  density  i-s  varyung  the  most. 

We  can  improve  on  this  estimate  as  foliow.s.  Suppose  in  the  neighbourhood  of  the 
maximum  density  gradient,  the  stratification  is  of  the  form 

i/(')  s;  -  piz  -  Zm'r 

for  some  positive  constant  p.  This  is  a  quite  general  relationship  since  i/(r)  must  have  a 
maximum  at  r  =  z^-  Then  suppose  -c  is  the  correction  to  i.e. 

k-  =  <7i/rr,  -  f 

Then,  writing  t  =  z  -  Zm-  the  above  ODE  becomes  approximately  (locally  to  t  —  0) 

I'E  -e  if  —  apt^  H'  =  0 

and  must  be  solved  subject  to  H'  — *  0  as  f  ^  ±5c. 

This  IS  a  standard  Hermite  or  parabolic  cylinder  function  ODE.  and  for  example  ha-; 
the  simple  exponential  solution 
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if  €'  =  afi.  More  generally  there  are  solutions,  which  are  polynoiuiaLs  of  degree  j  -  1  times 
the  above  exponential,  proriding 

r  =  (2 j  -  D'trp 

This  is  the  j'th  internal  wave. 

Thus  as  <7  —  oc . 

k~  -  -  {2j  -  Dv^  0(1) 

or 

k  =  v  -  U  -  '  0((T- ‘  - ) 

-  V  ^rn 

Hence  (asymptotically)  each  mode  is  obtained  from  the  previous  one  just  by  a  constant 
down  shift  in  k  of  magnitude 


This  answers  one  or  two  of  the  questions  in  Appendix  3.  In  particular,  in  the  last 
equation  of  that  appendix.  C  =  kj  —  Ak  (independent  of  mode  number  jj.  I 

don  t  know  what  is;  that  is  a  higher  order  effect  for  smooth  stratifications  such  that 
t^iz)  has  a  simple  maximum  at  an  interior  point  c  =  Experience  with  nonsmooth 
stratifications  (e.g.  stepwise  constant  i'(z))  or  with  those  where  the  maximum  density 
gradient  occurs  at  amend  point,  or  with  empirical  formulae  like  Philhps  s.  tends  to  suggest 
erroneously  that  an  apparent  rightshift  in  a  is  the  leading-order  correction,  rather  than  an 
apparent  downshift  in  k. 

[Figure  6.1  gives  samples  of  our  computation  of  the  k  -  a  curves. 
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Appendix  7:  11  March  1992. 


What  submarine? 

Most  of  our  work  so  far  has  ou  the  waves  made  by  a  siri';le'  isulatc'd  jxiiur  >eu;i  < 
of  uiiif  strength.  The  "fina]”  result  is  equation  (4-3)  of  Appendix  1,  iianudv 


ir  ~  -  /  sin{k.rcoi-^'H'o>{ki/>mH]A{kj(lH 

-  ./o 

where  Aik)  =  \Vc,k/Dk- 

At  this  stage  (Program  B)  we  are eoinputmg  .4(A-)  sun-ts:sful!y  for  tht>  unit  sourer,  fir 
various  stratifications  and  source  depths. 

Note  that  the  above  form  of  (4.3)  does  not  distinguish  niodrs  j.  the  results  brung 
assumed  computed  separately  for  each  mode  and  summed  at  the  end.  I  have  just  wTitten 
k  for  Kjia).  assuming  that  if  the  above  integral  is  being  earned  out  with  resw'ct  to  the 
theta  t'anable.  then  we  work  out  a  ~  xsec- tl.  then  call  upon  our  program  to  evaluate  the 
corresponding  k.  etc. 

A  serious  alternative  being  considered  at  the  moment  is  to  convert  directiy  from  a  V 
integral  to  a  k  integral,  i.e.  set  dB  =  IdB.’dk'idk.  working  out 


.2n  Si'c^  Bi>mB 
da 


which  demands  output  from  the  progiain  of  the  slope  lik  'da  of  ihv  dispersion  relation, 
which  is  not  too  hard  to  get. 

Importantly,  the  integral  with  respect  to  B  self- truncates  at  its  lower  end.  indeed 
only  occupies  a  small  range  near  B  =  r./2  for  realistic  small  stratifications.  That  is.  t.he 
generating  amplitude  A  is  identically  zero  for  a  <  a^  (the  lowest  raiue  of  a  for  each 
mode),  and  aj  is  a  large  number,  of  the  order  of  thousands.  Hence  the  actual  lower  limit 
of  integration  is  B  =  0,  where  =  k  set.-^  Bj  or 


which  is  very  close  to  ~,'2  unless  k  is  very  large  (slow  sub.  but  that  will  make  small  waves 
am-way) . 

This  is  all  well  known,  and  commented  upxm  by  mo.st  authors.  It  means  that  internal 
waves  are  largely  diverging  rather  than  tramsverse.  i.e,  their  crests  are  parallel  to  the  sub's 
track.  Bur  there  are  a  number  of  important  consequences. 

First  note  what  happens  when  we  make  the  change  of  variable  from  B  to  k.  Then 
the  range  of  k  is  fully  from  0  to  oc  (corresponding  to  B  =  B;  to  ~/2).  .Meanwhile  the 
X'Wi.se  wave  number  krosB  i.s  going  from  0  to  a  finite  upper  bound,  whose  magnitude  (see 
Appendix  6)  can  be  shown  to  be  a  small  quantity.  The  y-wise  wave  number  is 

essentially  k  itself,  and  goes  from  0  to  The  integral  converges  subject  to  only  mild 


re:r:trictious  on  .4  becau>t'  of  the  fact  that  d(i  dk-  tends  to  /.eio  Uke  k  ~  uvhu  ii  tdso  toilows 
from  Appendix  6  after  some  numipulatum). 

Consider  first  the  Raiikine  ovoid  reprr'scntation  of  a  real  suhimaine  ot  k  iiftth  L  and 
volume  r.  According  to  slender  boilv  thtxirv.  this  is  a  evlmder  of  uniform  cro»-st'cnou 
mea  V  L  =  irR-.  This  iiolds  so  long  as  R  «  L 

But  now  let  us  exanuue  what  hapjams  waen  we  start  to  evaluate  the  integral  for  the 
waves.  As  dist,'ussed  at  the  end  of  Appendix  5.  all  that  hap^KUis  relative  to  the  source 
integral  is  that  the  sin(A'.rcosi*^’  gef>  replaced  by  a  "cos  .  and  the  anijilitude  tai  tor  A  i- 
multiplied  by 

2sini  —  A•cu^ft! 


We  just  showed  that  the  largest  value  that  can  t>e  taken  by  t  he  .r-wist‘  wave  numbe:  k  eus  h 
is  the  small  quantit>'  Hence  the  argument  of  the  above  >.m  funcrioii  ran  ne\-er 

exceed  ^'Ku^Lr2. 

A  rough  estimate  for  realistic  st rat ihcatio us  i>  /y„  -  lb"  'm  '*.  >o  thn  (juantitv  romes 
TO  about  O.OOoi  T'  where  the  sutimarme  length  L  !>  in  metre>  mid  the  speed  ('  ui  n)etre> 
per  second.  1  believe  that  this  will  be  small  for  all  realistn-  -uiimarme  length.'  atni  <pf>t‘ds, 
at  least  small  enough  for  its  sine  to  Ih>  reasonably  replai-ed  b>'  it.'tdf. 

If  that  IS  done,  wtiat  results  is  the  same  as  if  wt'  had  let  L  —  b.  In  other  word.', 
all  submarmes  of  the  same  volume  generate  the  same  internal  waves,  irresjx’cnve  of  then 
length /diameter  ratio!  Hence  w  can  assume  the  submarine  is  a  sjihere.  winch  is  generated 
by  a  dipole. 

The  actual  effect  is  to  replace  the  sine  in  (4.d)  by  a  cosim'  and  t(j  multiply  the  lunph- 
fude  function  .4  by  the  .r-wise  wavenumber  A'cus6'.  and  the  final  resmr  bv  just  the  product 
of  the  speed  and  volume  of  the  submarine  UV.  Yet  anor  .  i  v.ac  'o  ,'.t  ,  :as  is  to  note  that 
the  waves  due  to  a  dipole  are  proportional  to  the  j  iciivanve  of  the  wavt's  due  to  a  soun  e. 

In  any  case,  we  should  compute  sot,  -  waves  due  to  dipoles  and  (separately  i  s>me  due 
to  Rankme  ovoids  of  the  same  volume,  and  set'  if  there  is  any  difference.  Dipohs.  hrst 


Appendix  8:  12  March  1992. 

A  review  of  rtie  DREP  report  by  T.\V.  Dawson 

The  1988  Canadian  rejxirt  DREP  TM8a-T  by  T  W.  Dawsoii  is  very  relevant  tu  (Jiir 
project.  This  note  is  a  sunuuar\-  of  some  unix)rraat  fear  tires  trf  it 

Overview' 

DREP  TM88-7  is  about  90  pages  long,  and  verv  detaile<i  It  an  impressive  docuiuenr 
in  its  accuracy  and  coverage.  It  reads  like  a  Ph.D  tiu-sis,  and  I  wouhl  rate  ir  ;is  (ieser\'ing 
of  a  Ph.D.  if  it  was  one.  On  that  basis  and  other^.  I  assess  its  cost  as  abont  .1  higli-kwe! 
man-years,  though  this  could  be  conservative  since  there  are  references  to  supporting  work 
that  was  obviously  done  in  parallel.  The  cost  of  the  rtsseaxch  dirts  tly  reported  in  TMxx-7 
alone  must  therefore  exceed  a  half  million  dollars;  if  the  supporting  work  is  als'  included, 
perhaps  the  total  cost  would  run  to  several  million  dollars.  Compare  our  fajiproximateiy  i 
6- week  times  4-person  effort' 

The  good  and  bad  news  is  that  this  repon  is  highly  .elevant  to  oiu  project  and  use>  a 
comparable  merhodolog}-.  I  do  not  believe  that  we  will  be  able  to  come  up  with  annlnng 
significantly  better  than  the  Canadians.  Hence  one  serious  alrernative  strateg%'  would 
have  been  simply  to  buy  in  their  finished  product.  If  that  was  not  iiossible.  we  could  have 
attempted  to  duplicate  it  using  TM88-7  as  our  guide,  I  hope'  to  make  it  clear  that  the 
latter  would  have  been  infeasible  in  the  time  available. 

That  is  probably  the  bad  news.  The  good  news  is  that  we  are  on  the  right  track  in 
our  independent  approach.  In  panicular.  I  am  absolutely  amazed  that  rny  prehminai-v 
note  •'Source  in  a  stratified  fluid"  (call  ir  EOT92-1  here:  it  is  reproduced  as  .Appendix  1 
of  the  final  report)  which  was  written  without  benefit  of  studying  TMt<8-7  hrst  describes 
a  methodology  that  is  quite  like  Dawson's.  There  are  many  other  things  oni'  could  do. 
and  some  in  our  group  thought  w'e  might  be  doing  them,  ranging  up  to  supercomputtT 
Xaviei-Stokes  solutions.  But  I  chose  a  more  conser\-ative  approach,  and  so  did  Dawson. 
We  could  both  be  wTong.  of  course,  but  there  is  some  experiment aJ  evidence  at  the  end  of 
IM  88-7  that  Dawson  s  results  are  good,  which  is  another  piece  of  good  news. 

Summary  of  TM88-7 

1.  “Introduction.’*  .A  nice  sumraaiy  of  the  problem.  Note  !pT)  the  discassion  of 
two  separate  generating  mechanisms,  "w'ake  collapse"  and  •'hull  effect"  and  an  important 
remark  that  the  methodology  allows  both  mechanism.s  to  be  studied,  although  the  illus- 
tration  is  only  for  the  second.  This  al,so  applies  to  our  approach.  The  approximations  and 
assumptions  (p.2)  are  as  I  w'ould  make,  and  (to  my  surprise,  since  I  thought  few  others 
had  even  heard  of  it!)  Dawson  also  chose  a  Rankine  ovoid  model  of  the  submarine 

2.  “Model  and  basic  equations.”  This  repeats  derivations  that  are  for  exainjiic 
given  in  Yih's  {196-5)  textbook  (.see  bibliography.  Appendix  2).  On  p.8  the  Boussinesq 
approximation  is  introduced.  I  do  not  propose  to  use  this  approximation  explicitly,  though 
I  am  sure  that  it  is  a  good  approximation  for  realistic  stratifications.  However,  as  a 
practical  matter,  if  does  not  provide  any  numerical  simplificarion  to  the  ODE  that  has  to 
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be  solved  for  the  vertical  distril)urion  of  veiociry.  and  hence  need  tiut  lie  made  excej>t  as 
an  analnic  tool.  The  comparison  is  betwet'n  iny  EOT92-1  etjuatiDii  (2  ')  i  and  Dawson  > 
TM88-7  equation  (3.4)  (subject  t-o  (3.2)).  My  txpiafion  reiluces  to  Dawson  s  if  '  ' 
is  replaced  b>-  pU'".  which  is  the  Boussinesq  approxiniatjon.  \iihd  in  the  iuun  of  samli 
stratification.  Note  that  his  .V-  =  -</p'  p.  and  iny  n  ~  </  r- 

3.  “General  solution  of  the  basic  equations.”  One  difft  rerue  our  ap¬ 

proaches  is  that  Dawson  uses  rigid  lid  conditions  H'  =  0  both  at  rlu'  -  ;;:.ii  e  :  --  (i 
and  the  "bottom"  :  =  d  (his  :  is  downward).  1  use  neither.  Mv  &•(>»’  .-urfacc'  is  gf-mnne  frev 
one.  where  W  =  a\V .  and  my  "bottom"  is  not  really  a  liottom.  Imt  rather  a  level  where 
the  density  becomes  uniform,  and  is  assumed  uniform  for  ever  below  in  an  inhaite  depth 
ocean.  Xeither  of  these  differences  is  too  imjxirtaut  m  practice  The  a[)parent  ngiduv  of 
the  free  surface  is  consistent  with  the  Boussinesq  approximation:  the  liortom  condition  is 
probably  not  too  important  in  practice  unless  the  sulnnarme  is  close  to  i? 

Dawson's  equation  (3.11)  looks  familiar  to  me.  It  is  closely  connix-ted  witli  mv  EUTd2- 
1  (4.3).  though  arising  very  differently.  In  both  <-ases.  the  Sfrnmia  involves  m  us  d^"!iomlna- 
tor.  the  partial  derivative  with  respect  to  a  wave  pariuneter  of  tlu‘  WroiLskian  Ix-twemi  two 
solutioms  of  the  fundamental  homogeneous  ODE.  Tht'se  solutions  are  m  geiu'ral  linearly 
independent  (so  the  Wroaskian  is  non-zero)  but  are  of  mam  interest  to  us  when  depeiulenr 
(so  the  Wronskian  is  zero).  In  that  case,  the  wave  number  k  and  mv  parameter  n  leepiiv- 
aiently  the  wave  speed  c)  are  connected  numerically  by  tht'  c  spersiou  relation.  Dawsoi, 
avoids  actual  differentiation  of  the  Wronskian  by  some  mampulation  momul  p.  I'J.  but  ! 
don't  see  any  problem  doing  it  by  brute  force  numerical  differentiation. 

The  summation  over  modes  in  (3.11)  and  throughout  TMf'3-7  represeuits  the  major 
point  of  departure  between  the  philosophies  of  Dawson  and  me  (and  1  think  a  number 
of  others  are  with  me  on  this).  TM8k-7  adopts  the  view  that  one  must  compute  vuiny 
internal  waves.  There  are  infinitely  many.  They  are  ordered  with  resjieet  to  s{xvd  riare 
being  a  maximum-speed  mode  j  =  1  (minimum  a  in  mv  notation),  then  a  slower  mode 
j  =  2  (larger  a)  etc.  Roughly  speaking,  the  speeds  go  down  like  the  inver'^e  of  the  moije 
number  j  (ct  goes  up  like  f-).  Since  a  is  already  rather  large  even  for  the  first  mode  j  -  1. 
the  actual  <t  numbers  for  higher  modes  pretty  quickly  get  huge.  It  is  rm-  view  that  almost 
all  that  one  wants  to  know  about  internal  wave  generation  is  available  by  assuming  that 
only  the  first  mode  is  generated.  We  can  calculate  several  modes  if  we  want  to:  rh(*  results 
are  in  any  case  just  additive,  so  we  can  test  this  hypothesis  easily:  I  wonder  if  Dawson  did'.' 

4.  “Uniformly  translating  source:  frequency  inversion.”  Dawson  ba-  a'-surnt'd 
a  general  source  distribution  doing  somewhat  general  things  up  to  now.  Now  he  speriali.>es 
to  -Steady  horizontal  tramsiation,  as  I  have  done  from  the  out.set.  mid  first  does  an  isolated 
source. 

5.  “Extraction  of  steady-state  fields.”  More  specialisation. 

6.  “Evaluation  of  integrals”.  It  appears  that  Dawson  is  mnbitious  enough  to 
try  to  evaluate  the  whole  thing.  Thar  is.  in  my  terminologv-.  he  ecaiuates  not  only  the 
free  wave  single-integral  contribution  from  the  residues  at  the  poles  k  -  Kj{/r).  hut  also 
the  doibie  integral  local  disrurbancf>s  near  the  source.  .As  I  say  on  p.  4  of  EOTr)2-l.  'hi- 
IS  a  truly  ■niarmnoth  task".  I  can  only  have  admiration  for  Dawson  for  attempting  it 


aiid  more.  Perhaps  we  in  the  shi{>-hydrodynai«K-  community  iiave  b(>t'n  too  coasetAHtive, 
but  I  do  have  my  doubts  about  whether  one  can  really  do  it  with  accuracy  and  etiicKUicv 
J.N.  Xewman  has  been  quite  scathing  in  the  past  about  some  attempts.  Aiivwav  if  one 
is  mainly  interested  in  the  far  field,  the  free- wave  integriil  EOT92-1  should  suffice. 
Here  I  am  calling  upon  decades  of  ship- hydrodynamic  experience,  with  a  slight  wony  that 
some  of  this  experience  might  not  translate  dirtvfly  to  internal  wa\es. 

7.  "Ntimerical  consideratjons”.  This  is  a  very  long  and  detailed  chapter.  It 
describes  Dawson's  technique  for  solving  the  ODE.  his  (3.4).  ('qinvalent  to  EOT92-1 
N'amely.  simply  represent  the  ocean  by  many  layers,  in  each  of  which  the  density  is  assumed 
to  vaiy  exponentially,  so  A'*  is  constant.  Fit  thtsie  layers  together  with  suitable  continintx 
conditions,  and  you  have  a  tool  for  ODE  solution. 

Much  of  this  chapter  is  taken  up  with  concern  for  the  accuracy  of  this  techiuqm 
related  to  things  like  "evTinescent  modes"  (complex  eigenvalues),  mid  close  ajiproach  of  the 
real  eigenvalues  for  one  mode  to  those  for  another. 

I  believe  to  a  certain  extent  these  worries  are  an  artefact  of  the  nif'thod  Dawson  is 
using  to  solve  the  ODE.  though  they  are  a  warning  to  ns  also.  We  are  using  a  direct 
nuraerical  solution  method  of  Runge-Kutta  type  (not  for  any  rt'ason  other  than  ea-e  of 
programming)  and  I  don't  think  we  are  seeing  any  similar  difficulties 

Figure  2  on  page  38  is  very  important.  It  provides  a  benchmark  density  distribution 
(indirectly  via  the  Brunt-Vaisala  frequency  .V  =  y'^).  We  are  using  it  to  rix'over  the  actual 
density  p  by  integration  and  exponentiation,  and  then  inputting  that  to  our  program.  Thi> 
is  a  bit  silly,  since  in  effect  our  program  just  differentiates  p  immediately  to  give  n  again, 
but  is  justifiable  on  convenience  grounds.  After  some  effort.  I  think  we  have  quite  good 
data  for  p  at  intervals  of  10  metres,  as  used  by  Dawson  in  his  Figure  2.  and  also  interpolated 
to  5  metre  intert'ais.  This  density  data  is  shown  as  a  Figurt'  here 

Dawson's  Figure  3  on  p.41  is  one  that  we  should  be  able  to  reproduce.  It  is  m  effect 
the  dispersion  relation  plot,  equivalent  to  our  k  versus  tr  plots.  IndetxJ  there  is  considerable 
equivalence.  Dawson  s  variable  s  is  proportional  to  our  n.  and  (set'  his  e<p  I.13),i  hi.';  p 
is  proportional  to  :r;.  So.  apart  from  stjrne  possible  confusion  viver  the  proportionality 
coefficients,  once  we  find  k  versus  a.  we  have  his  Figure  3. 

Of  course,  there  are  21  modes  on  his  Figure  3;  we  are  rather  more  modestly  thinking 
of  computing  about  3  modes!  One  difficulty  I  find  with  Figure  3  is  that  1  can  r  tell  which 
mode  is  which.  The  lowest  curves  should  be  the  earliest  modes,  but  the  first  3  of  these 
are  incomplete.  The  logarithmic  scale  doesn  t  help.  Am-v^ay.  the  sooner  we  attempt  this 
comparison  the  better. 

The  rest  of  Chapter  7  is  special  nuraerical  techniques,  and  perhaps  not  too  iriteresnng. 
Figure  7  (p.64)  starts  to  be  interesting  again,  and  all  the  remaining  Figure's  should  be 
reproducabie  by  u.s. 

8.  “Wake  examples”.  Here  on  p,71  Dawson  introduces  the  Rankme  ovoid  idea, 
and  gives  a  concrete  example  with  submarme-like  dimensions.  He  then  gives  final  rcsulrs 
for  that  example,  at  several  speeds  and  depths  of  submersion,  in  the  stratification  of  his 
Figure  2.  Figures  12-14  are  3D  plots  which  we  should  be  able  to  reproduce  of  the  surface 
currents.  However,  the  actual  scale  of  these  plots  is  not  easy  to  jiick  up.  though  there  are 
Tables  of  extreme  values  provided. 


In  the  caption  to  Figure  12.  the  ‘iocal"  effect  near  the  submarine  is  noted  We  can't 
reproduce  that,  since  we  neglect  the  local  effect  via  the  double  integral  term  m  the  wave 
field.  If  the  local  effect  is  as  small  as  that  shown  in  Figure  12.  this  'liould  not  be  too  much 
of  a  woria-. 

9.  “Comparison  with  measured  data”.  Figure  Iti  is  a  very  mtens^nng  compm- 
ison  between  computation  and  experiniem.  the  experimental  internaj  wavis-  being  shije 
generated.  The  agreement  is  very  encouraging,  both  qualitafi%T'ly  and  quiUitiranvtUw  We 
could  not  directly  reproduce  this  without  nutre  informanoti.  presumalily  obtainable  from 
Reference  1.  which  is  a  DREP  report  of  19tf5.  I  would  like'  to  know  how  nuui\’  internai 
wave  modes  were  significant,  and  whether  "lot'a!''  effects  dlK'  double  integral  terms;  wire 
significant. 

10.4-  “Concluding  remarks.  Reference  and  Distribution  Lists."  I  '.s<-ful  The 
US  reports  under  the  name  Milder  seem  likely  to  l)e  imertssting.  There  are  some  interesting 
inclusions  (Cart'RTighr.  Reed)  and  omi.ssions  (Tulin)  in  the  distribution  list. 
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Appendix  9:  13  March  1992 


Direct  ^^-integration 

After  some  thought.  I  have  decided  that  (at  lea^t  for  tJie  time  being)  we  will  du  ihe 
$  integration  directly,  i.e.  without  the  change  of  variable  to  k  ranvassi’d  m  Appendix  7 
The  alternative  is  still  worth  considering,  but  let  s  go  with  t^-imegration  first. 

Hence  the  vertical  velocity  made  by  a  Rankme  ovoid  of  length  L  and  iiuLxjinmn  radiu.- 
R  moving  at  speed  V  is  given  b>‘ 

2V ’  f~  '  siiH  k  (■( IS  fi'i  -  2 !  1  ,  . 

- /  cos(A'xcost^icos(A-vsintt,i - ^ - AikUift 

where  V  =  t:R~L  is  its  volume. 

The  lower  limit  of  integration  vs  Sm-  which  is  the  nhiumuni  «  value,  corresponding  to 
the  minimum  a  —  Um  value  for  each  (separate)  internal  wave  mode  j  Xaniely 


um  —  ^ 

\ 

We  assume  for  validity  of  this  formula  that  k  <  am-  which  is  niit  a  severe  requirement  since 
(Tm  is  numerically  large.  This  is  what  Tulin  and  Miloh  (1990:  se<>  bibliography.  Appeiuiix  2) 
call  (misleadingly)  the  "supersonic"  range.  This  just  means  that  the  submarine  is  moving 
faster  than  the  fastest  internal  wave,  namely  faster  than  aln^ut  one  metre  per  si'cond. 
Absence  of  smaiU^  contributions  means  absence  of  transverse  waves.  If  the  sufimarine  is 
moving  so  slowly  that  this  condition  is  violated  ("subsonic"  motion  in  the  sen,s«‘  of  Tulin 
and  Miloh  1990).  the  generated  wave  amplitude  is  probably  too  low  to  detect  aip'way. 
However,  if  we  did  need  to  compute  it.  we  would  just  set  -■  0.  since  then  all  vaiiit's 
occur,  and  there  are  transverse  waves. 

Since  is  large.  Bm  is  very  close  to  r/2.  The  range  of  integration  abo\'e  is  therefore 
a  very  narrow  one  with  respect  to  9.  I  think  this  is  good  news  for  convergence  of  the 
integration.  Actually,  that  convergence  is  also  illustrated  by  the  k  transformation,  which 
leads  to  an  integration  on  an  infinite  A-range.  but  with  A-'*  rare  of  decay  of  the  integrand. 
In  some  sense  this  decay  rate  with  respect  to  k  is  equivalent  to  the  small  w-range. 

Traditionally  (for  surface  waves)  we  are  worried  in  ev-aluaring  this  son  of  integral 
about  the  effect  on  our  integration  method  of  a  rapidly- varying  character  to  the  integrand. 
This  comes  from  the  trigonometric  functions  when  x  and  y  are  large,  strictly  when  kx  and 
ky  are  large.  That  occurs  when  we  are  in  the  far  field  of  the  disturbance,  strictly  many 
wavelengths  from  it. 

There  are  several  reasoas  to  suspect  that  this  won't  be  a  problem  in  oiu  case.  First .  the 
internal  wavelengths  are  much  greater  than  those  for  surface  weaves,  so  at  a  fixed  (p>erhaps 
large)  x.  y.  the  value  of  kx.  ky  may  not  be  so  large.  Second,  there  is  an  actual  upper  bound 
on  the  coefficient  of  x.  namely  krosB  <  (see  Appendix  7)  which  is  small.  Hence 

there  cannot  be  a  rapidly  varying  character  to  the  contribution  from  the  factor  m  x.  no 
matter  how  large  x  is  Since  9  is  close  to  t,/2.  the  y  factor  is  es.spntially  cosily),  and  will  be 


rapidly  varying  when  y  is;  large.  However,  we  are  rather  less  inten^sreti  lu  litrgf  y  hti  tut  h*- 
side  of  the  submarinel  than  we  are  ui  large  .r  (far  ftstern  of  the  >!ihmarutei  Fiuaiis  !  iu-rt  s- 
only  one  point  of  stationary  pha-;^’  (so  long  a-  ^  )  and  Miat  i'  eioM  r,  i  -  j  .a  ri 

since  all  B  values  have  that  property.)  If  we  were  to  be  vvornui  about  the  iioigiuM  .urhoo.i 
of  points  of  stationary  phase’,  we  would  Iw’  worrvmg  about  rlu*  luugliiHHn  iiood  ot  e  - 
where  there  is  no  contribution  at  all  from  the  .r  factor  Tlie  la.":  .iigtimenr  i-  [notiabh 
facile;  but  an\-way,  all  1  am  arguing  is  that  we  should  pres.'  ahead  doing  riif  integr.ition 
without  worrying  too  much,  and  the  proof  of  the'  pudding  i^  in  tlie  eating. 

In  Appendix  7.  it  was  argued  that,  as  gcnerattirs  of  internai  wave>  all  Miiun.trini''  ot 
the  same  volume  are  indistinguishable.  irres{)ect iw  of  sha{H’  mid  length  I  tluuk  tin.'  S'  .m 
important  principle  that  we  can  check,  nuxst  directly  for  the  Raiikine  bodK"  b\  varving 
their  length/diameter  ratio  L  {2R)-  In  particular,  in  the  hnur  as  hL  --  0  ‘note  length 
tending  to  zero  relative  to  wavelength,  nor  necessariiy  relative  to  diameter/  we  itsover 
the  internal  waves  generated  by  a  sphert'.  where  the  quanntv  in  >quare  bracket>  is  'impiv 
replaced  by  A'ccsW.  and  the  result  is  independent  of  L.  It  would  be  very  interesting  to 
check  this  by  computations  at  a  sequence  of  L  values 

Some  comment  on  dimensions  is  appropriate  here.  Obvioiisiv  the  output  a-  C  must 
he  dimensionless.  The  volume  factor  outside  the  integral  has  dimensions  uf  length  l  ubed. 
so  the  integrand  must  have  dimensions  of  length  to  the  pt)wer  -  d.  Tin.'  is  true  since  n  is 
not  hard  to  see  that  .4  has  dimensions  of  length  to  the  power  -2  (.se<*  .^ppa-ndix  1.  below 
(3.4);  Uo/D  is  dimensionless*. 

I  suggest  that  we  use  consistent  SI  dimernsions.  i.e.  input  and  use  L  in  metres  etc 
Up  till  now  we  have  normalised  all  variables  so  that  the  apparent  deptli  of  the  ocean  H  i.' 
unity.  I  think  this  is  no  longer  appropriate  in  the  production  program,  .so  we  should  input 
the  actual  depth  H  in  metres,  etc.  Since  all  results  are  proportional  to  the  submarine 
volume  V,  that  factor  need  not  be  inputted.  If  it  is  omitted,  we  are  computing  the  wave 
made  by  a  submarine  of  volume  one  cubic  metre;  a  submarine  of  volume  l(M)b  cubic  metres 
makes  KKK)  times  that  wave.  Otherwise,  all  variables  are  in  true  SI  units. 
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Appendix  10:  17  March  1992 

Density  data 

As  seen  in  Appendix  GF  (a  preliminan'  note  by  Graham  Furnell).  tiie  data  is  not 
always  presented  to  us  directly  in  the  form  of  density  versus  depth  pi:].  It  may  lie  salmitv- 
or  temperature  profiles,  or  it  may  be  speed  of  sound  data  or  even  Brunt-\'aisala  freipiencies 
whose  square  is  proportional  to  our  intermediate  variable 

Generally  though,  we  expect  to  see  a  p{:)  that  increases  with  depth  (aimliing  else 
being  statically  unstable)  and  our  results  only  depend  on  the  ratio  between  two  densities 
so  without  loss  of  generality  we  can  assume  that  the  densit\  at  the  free  surface  i.-  normalised 
to  unity.  Then  it  will  increase  with  depth  and  the  actual  mcrt-ases  are  generallv  less  riuui 
one  percent. 

The  action  takes  place  in  the  first  few  hundreds  of  metres.  There  may  be  a  mixing  zone 
for  the  first  few  tens  of  metres  where  p  is  essentially  constant,  then  {auhaps  a  relatively 
rapid  increase  with  a  peak  density  gradient  in  the  next  few  teiLs  of  metric  {the  t  hermoclme  i . 
then  less  dramatic  increases  in  density,  with  a  levelling  off  at  about  metrts^  to  a 

final  density  of  =  l.OOlo  to  1.0070,  (times  the  fret'  surface  value)  at  least  in  the  examples 
I  have  seen 

There  are  interesting  examples  where  the  Brunt-Vaisala  frequency  .%'(:*  is  quottd  m 
radians  per  second.  [Sometimes  it  is  given  in  cvcies  fier  htiur.  which  must  fie  first  convert txi 
by  multiplying  by  2r/3600|.  One  way  to  proceed  with  these  would  be  to  put  this  data 
directly  into  our  program.  Instead,  however,  it  is  (lossibie  to  recover  tht'  actual  dert'.itv 
profile  from  .V(-)- 

Namely,  since 

=  -gp’-p  =  9<' 

then  (normalised  to  unity  at  c  =  0) 


pi:} 


X'iti/gdt) 


We  have  done  this  integration  for  a  few  .sample  data  sets  The  most  iitifxirtant  is 
Dawson's  Figure  2.  reproduced  here  as  Figure  10.1.  The  consequtmt  densitv  curve  is 
Figure  10.2.  A  decision  has  to  be  made  as  to  where  to  cut  off  the  data.  Dawson  shows  .V 
as  non-zero  even  down  to  depths  of  the  order  of  kilometres.  It  is  a  good  tiuestion  whether 
this  data  is  belie\able.  but  in  any  case,  the  actual  values  of  .\^  are  small,  of  the  order  of 
hundreds  of  times  less  than  the  maximum  of  so  the  density  is  only  increasing  very 
slowly  at  great  depth.  We  have  to  do  some  numerical  experimentation  to  verifr  this,  but 
I  would  expect  that  we  would  be  making  very  little  error  if  we  were  to  assume  that  the 
density  was  constant  w'henever  2  is  greater  than  about  300  metre's  in  the  conditions  of 
Dawson's  data,  and  also  for  .some  other  examples  we  have  seen.  There  are  some  cases 
where  there  does  seem  to  fx'  significant  variation  down  to  order  of  a  kilometre,  but  I  doubt 
its  significance  for  our  problem. 
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Stobiiity  Profile 


3la  frequency  profile  used  for  the  examples 
smoothed  version  of  measured  data  taken  at 


mm 


Appendix  11:  20  March  1992 

Early  output 

We  are  now  in  a  position  to  produce  some  preliminarv'  i)utput  q\umt  it  tliat  can  and 
in  some  cases  have  been  checked  against  existing  results. 

In  the  first  place,  regarding  just  the  kinematics,  namely  the  di.>:persion  relation,  we 
have  checked  a  number  of  Dawson's  results,  using  our  own  n'-coast ruction  (Figure  ib.2)  of 
his  density  profile.  That  reconstruction  is  a  potential  source  of  error,  but  the  result.'  s«>r« 
pretty  good,  so  I  think  our  reconstruction  must  be  OK. 

For  exanrple.  Figure  11.1  reproduces  Dawson's  Figure  'i.  as  anticipattHl  in  .^ppeiuii.x 
8.  We  have  to  pick  a  value  for  "-Vn,ax  •  ^  TM88-7.  [).  d7,  .According  to  Figure 

2  of  TM88-7,  A'sna.x  =  14.6  cycles  per  hour  •,  at  least  that  is  the  largest  c-alue  of  the 
Brunt'Vaisala  freo.uency  displayed  on  that  .step-like  Figure,  and  we  believe  it  is  likely  that 
that  was  the  value  used  by  Dawson  for  his  calculations  leading  to  Figure  d.  Interestingly, 
however,  our  first  choice  was  to  '’smooth"  the  steps  in  Figure  2.  which  leads  to  a  "better " 
estimate  A’max  =  15-9  cycles  per  hour,  since  there  is  a  very  sharp  maxtmum  that  is  nut 
well  captured  by  the  steps  in  Figure  2.  The  result  was  tt)  shift  all  our  results  by  a  small 
constant  amount,  which  (on  the  logarithmic  scale  used)  corresponded  to  the  approximately 
207c  factor  in  A'^^.x- 

In  any  case,  with  what  we  beheve  was  Dawson's  choice  of  Am^x  axe  always  spot  r)n 
his  graph  for  modes  1  and  2.  Dawson  gives  21  modes,  but  only  incomplete  graphs  for  the 
most  important  first  3  of  t'.ese.  We  .show  (dashed)  the  complete  graphs  for  modes  1.2  and 
3:  those  for  mode^  j  2  merge  smoothly  with  Dawson's  and  extend  Dawson's  rf^sults 
correctly  toward  p  i.e.  toward  k  ~0 

The  sam  me  for  mode  3  with  the  results  shown  in  Figure  11.1.  which  wa<  bas^i 
on  use  of  d^'nsity  data  extending  down  to  about  a  kilometre  depth.  However,  if  is  worth 
noting  That  our  first  computations  using  data  extending  down  to  oniv  500  metre's  depth 
diverged  from  Dawson's  mode-3  cmve  near  its  bottom  end.  This  is  because  mode  3  has 
a  maximum  at  about  600-800  metres,  and  is  not  adequately  captured  by  500  metre  data, 
^lowever.  the  actual  generated  amplitude  (see  the  sample  curve  of  Figure  11.7)  is  smaller 
by  a  factor  of  about  10  than  the  corresponding  mode  I  and  2  signal  reflecting  the  fnct 
that  a  submarine  at  about  100  metres  depth  is  not  an  efficient  exi’iter  of  a  mode  whose 
maximum  is  much  deeper. 

As  of  this  date,  we  are  rapidly  producing  graphs  of  A{z.9).  Tht'se  give  a  very  good 
picture  of  dependence  on  parameters  like  submarine  depth,  mode  number,  for  a  parricuiar 
stratification.  We  are  presently  using  Dawson's,  but  that  is  a  bit  weird,  and  wc'  should  shift 
to  some  more  realistic  stratifications.  The  Figures  to  follow  axe  prehminary  indications  of 
behaviour  oi  A(z)  for  various  a  (equivalently  k)  values.  The  first  4  Figures  show  mode-1, 
for  various  submarine  depths.  Figure  11.6  is  a  mode-2  for  depth  100  metres,  and  Figure 
11.7  mode-3  for  depth  100  metres, 
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A(  K(Sigma),  Z  )  for  the  1st  Internal  Wave 

Profile  =  dawson260.pro  Submarine  Depth  =  150.00 


Figure  114:  Mode-1  ar  h  1.^0 
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Figure  11,5:  xMode-1  at  h  =  200 
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A(K(Sigma),Z)  for  3rd  Internal  wave 

Profile  =  dawson  lOOO.prolO  Subm-Tine  Depth  =  280.00 
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Appendix  12:  23  March  1992. 

Other  components  of  velocity,  and  large  wavelength. 

For  reasons  associated  with  the  e.xact  linear  theory  (st'e  .\ppendix  I  )  we  have  con¬ 
centrated  on  the  vertical  component  «•  so  far.  A  somewhat  disiippomtiu<!.  a-]>cc'  of  tiiat 
concentration  is  that  the  \'alue  of  this  component  is  very  small  at  the  fr<x'  siirr'ace  I?  i- 
not  actually  zero,  but  appears  so  on  the  plots,  relative  ti*  rhf'  value>  at  dc'pth?-  o!  i  he  order 
of  100  metres  or  so. 

In  fact,  it  is  consistent  with  the  Boussinescj  approximation  that  u-  ^  0  at  th"  friv 
surface,  a  result  that  can  be  seen  by  noting  the  linear  free  surfaci-  condition  U"  ~  oU' m 
the  limit  of  large  a. 

From  the  practical  point  of  view,  the  other  two  vek)Cit\-  componenr.'  u.  '■  me  iiujre 
significant,  and  their  \'alues  are  at  least  as  groat  in  order  of  masmtude  at  th.e  fr«x'  siirface 
as  they  are  at  depth.  These  two  components  are  the  gradient  of  a  potent iai  with  Fourier 
transform  4>  available  from  the  equation  above  (2.3)  of  .\ppendix  1.  namely  4>  - 
Hence  as  well  as  the  leading  equation  of  Appendix  9.  namel\-. 
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wnere  B  =  A' ;k  =  W'^/Dk-  the  prune  as  usual  meaning  d  Oz. 

That  is.  to  compute  these  velocity  components.  ju.st  substiture  the  ;  dcnvanve  U'f  for 
the  function  Uq.  and  delete  the  factor  k  in  the  original  definition  .4  =  A-Ui,  Dc  to  give  a 
new  amplitude  function  B.  Then  replace  .4  by  B  in  the  integrals,  and  for  u  iucludf'  a  fartor 
-(•os6  and  for  v  a  factor  -sm6.  Finally,  the  r  and  y  factors  are  swapfxxi  appropriately 
between  cosines  and  sines.  The  result  is  no  more  difficiilt  to  compute  than  tr. 

Estimates  of  orders  of  magnitude  are  interesting.  First  note  that  u  is  the  disturbance 
T-component  of  velocity  due  to  the  submarine,  the  total  being  U  ^  u  m  a  frame  of  referfuice 
moving  with  the  submarine.  Then,  at  least  in  the  far  ■sufx^rcritica]"  speed  range  when 
K  <<  f7j  and  hence  8m  is  close  to  77/2.  we  must  have  w  <<  r.  since  tlien  cu>8  «, 
sinf^.  Thus  the  lateral  velocity  disturbance  should  he  more  significant  than  the  along- 
track  disttirbance.  This  does  seem  to  be  m  rough  agreement  with  Dawson's  resultx. 

.4t  the  same  time,  there  are  a  couple  of  arguments  .mggestirig  that  r  and  w  are  of 
the  same  order  of  magnitude.  Recalling  that  w  is  essentially  zero  at  tiie  fret-  surfact' 
whereas  i-  is  not.  this  means  that  the  value  of  r  at  the  free  surface  may  be  comparable  in 
magnitude  with  that  of  \i:  at  depths  like  100  metres.  In  the  first  place,  we  may  consider 
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that  the  important  values  of  the  wavelenRth  2r:  'k  iire  e-iiniparalhe  with  sigmticaiit  liepth 
scales  for  rates  of  change  of  the  density  stratification.  Hence  the  order  of  iiiagiutude  of  k 
IS  comparable  to  that  for  d  dz  Hence  B  and  A  should  have  conip;u'ai)le  rnagiutude>;  aiul 
thus  so  do  c  and  w. 

Another  related  intuitive  idea  is  expressexl  by  the  attached  Figure  12.1  which  is  not 
new.  but  appears  in  other  publicafioiLs  on  internal  wavt-s  generated  bv  steaiiiiy  moving 
lx)dies.  Namely,  when  we  have  neglected  u.  continuity  in  the  cross-how  juiuie  deuiaiids 
that  a  snapshot  of  the  flow  field  takes  tiie  form  of  clo.sed  loojis 

I  do  not  believe  that  this  means  that  Hind  particks  nei-issm-ily  undergo  such  loops, 
there  would  not  in  any  case  be  time  for  particles  to  move  round  such  loop>  the  sutunaruie 
passed.  This  matter  calls  for  further  study.  But  nevertheless  one  get?  the  impression  from 
such  pictures  that  there  must  lie  a  rough  balance  lietwwn  the  magnitudes  of  r  ami  n 

Wave- like  behaviour  occurs  in  each  of  the  components  r.  a-  as  functions  of  v  and  these 
waves  are  hkely  to  be  nearly  out  of  phase.  Hence  when  w  is  at  it?  greatt'sr  upward  value  at 
depth,  meanwhile  r  should  be  nearly  zero.  In  Itetwtxui  two  such  values  of  //  there  will  be  a 
near-zero  point  for  w  (at  depth:  ir  is  always  near  zero  at  the  surface).  At  that  laiue  of 
y.  then  r  will  be  near  maximum,  and  that  will  be  esixxially  true  at  the  frtx’  -u  face  i'si'if. 
etc. 

Figure  12.1  apphes  to  the  first  internal  wave  nuxle  only,  and  the  loojis  extend  down 
to  the  effective  end  of  stratification  for  that  mode.  A  similar  picture  mu.st  liold  for  higher 
modes,  e.g.  Figure  12.2  for  mode  2.  with  two  liKips  in  the  depth  dire<-tion.  The'C  loops 
are  pure  speculation.  We  can  compute  them,  though,  and  should  do  s<.i  ASAP 

It  is  always  well  to  keep  in  mind  that  the  general  length  scale  of  interest  for  tlu'se 
internal  waves  is  that  of  the  stratification,  namely  the  order  of  about  one  hundred  metres. 
This  scale  starts  life  as  a  scale  in  the  r-direction.  and  always  correctly  indicates  that  scale, 
so  that  any  disturbance  penetrates  up  and  down  over  these  rather  largt-  disTanco  with 
ease.  In  particular,  whatever  internal  wave  is  present  at  the  surface  will  al>o  !>e  presmit 
at  several  metres  depth  with  essentially  no  change  in  amplitude,  and  will  only  actually 
attenuate  at  depths  of  the  order  of  many  hundreds  of  metre's. 

Meanwhile,  this  scale  is  also  determining  the  horizontal  length  scale-  of  the  wave  Other 
factors  (such  as  the  u.  r  disparity  above)  suggest  some  disparity  lietwe'en  it?  manifestation 
in  the  X  and  y  direction-S,  Roughly  speaking.  I  fieiieve  that  the  y-wise  wavf'length  prese-n-es 
the  depth  scale  of  about  a  hundred  metres.  But  the  .r-scale  may  be  stretcheeJ  to  sf'veral 
hundreds  of  metres,  so  that  the  internal  wave  appears  to  be  propagating  mainly  sideways 
to  the  submarine's  track:  a  mainly  "diverging"  wave  pattern.  In  any  case,  these  are  very 
long  waves,  compared  for  example  to  surface  waves  made  by  a  ship  for  even  a  submarine), 
which  have  wavelengths  of  the  order  of  only  several  metres,  nearly  two  orders  of  magniTude 
smaller. 

Internal  waves  are  thus  qualitatively  unlike  ship  waves:  they  could  be  said  to  look  more 
like  "currents'.  They  are  unlikely  to  be  observed  directly,  but  rather  via  their  influence  oii 
other  observable  quantities  like  capillarv  wa\-es.  Indeed,  one  (perhaps  controversial'.')  sug¬ 
gestion  has  been  that  these  roughness-like  elements  of  the  micro-scale  frfx-  surface  pat tf'rn 
become  bunched  together  (higher?)  where  there  is  an  up-current  m  Figure  12.1.  i.e  where 
w  has  a  (yiositive)  maximum,  and  bwome  lower  in  inagnitudi'  where  there  is  a  (negative) 
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Further  effort 

1.  Mode  number  effects. 

I  would  have  liked  to  do  more  comparative  wiirk  on  the  inHuence  of  the  various  internal 
wave  inodes.  Dtx>s  the  generating  efficienc\-  dtH'rease  sigmficant  1\-  with  niotle  nunilier'.' 
How  many  modes  does  one  in  practice  have  to  consider'.’  How  dttes  the  nunifuical  acciir<u'\’ 
vary  with  mode  number  -  since  there  is  more  strucnue  with  n>i)ect  to  depth  m  the  higher 
modes,  do  we  need  a  finer  numerical  grid  for  them'.'  The  present  program  is  rather  primitive 
in  the  way  that  it  picJcs  out  the  particular  mode  that  is  desirtnl.  It  could  be  unproved  in 
that  respect. 

2.  Fourier  outputs. 

The  quantity  A{k)  on  which  we  have  spent  most  time  is  a  measure  of  wave-generating 
efficiency  in  the  space  of  the  wave-number  k.  However,  tuir  use  of  it  has  b^ni  directed 
toward  subsequent  integration  to  give  the  actual  point-wise  flow  field  iriJ.  y.  - ).  There  are 
other  possibilities.  In  panicular.  there  may  be  ways  to  detect  A{k)  ir.self.  m  t'ffect.  Some 
of  these  ideas  are  contained  in  the  Tuck.  Collins.  Wells  (1971)  paper  (sw  .\pf)eiidix  2). 
and  essentially  I  am  suggesting  that  the  present  work  be  combined  with  that 

3.  Kelvin  angles 

.■\nother  related  matter  that  we  gave  some  thought  to  early  in  the  project  but  neve'r 
found  the  time  to  pursue  further,  is  the  simple  Kelvin  stationary  phase  angle.  This  is 
discmssed  in  .Appendix  1.  and  there  is  no  reason  why  the  *  versus  6*  curves  could  not  be 
computed  routinely  as  additional  output  from  our  program.  There  is  a  separate  ciirvi'  for 
each  internal  wave  mode  and  (for  internal  waves  as  distinct  from  surface  waves  in  water 
of  infinite  depth)  a  distinct  curve  for  every  separate  speed  of  the  disturber.  1  believe  that 
the  results  may  be  in  some  way  similar  to  those  for  the  Kelvin  waves  in  shallow  water,  at 
supercritical  speeds,  for  which  (as  part  of  my  seminar)  I  drew  up  the  attai.-hed  curve. 

4.  Other  stratifications 

We  have  given  perhaps  unwarranted  attention  to  a  particular  stratification,  namely 
that  used  by  Daw-son.  The  only  reason  for  this  is  so  that  we  can  compare  wuth  some  of 
Dawson's  results.  But  that  stratification  is  nor  typiccii.  and  is  perhaps  somew'hat  bizarre. 
There  are  many  other  more  typical  stratifications  that  should  be  tried,  including  some  of 
direct  interest  for  the  present  application  to  detection.  A  related  numerical  matter  is  the 
question  of  how  deep  to  truncate  the  data.  We  have  used  data  from  Dawson  of  up  to  4t<0 
metres,  but  this  seems  urmecessary:  only  the  first  300  metres  is  really  significant,  but  more 
testing  is  needed. 

5.  Dispersion  relations 

.Appendices  3  and  6  contain  ideas  for  analvtic  work  of  a  novel  nature  on  the  character 
of  the  dispersion  relations  for  internal  waves,  .Just  what  features  of  the  stratification 
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do  the  parameters  of  the  dispersion  relation  depend  upon'  Pliilhp>V  formula  Mity;e>Ts 
things  about  the  separate  dependeuee  on  depth  and  thickness  of  the'  thermoclme.  iuit 
this  depends  on  a  particular  assumption  about  the  stratification,  and  riie  whole  idea  of 
a  "depth"  or  a  "thickness"  or  even  a  •Thennocline"  is  not  necessarilv  nieaiimgful  for  a 
general  stratification.  The  ideas  in  these  Appendices  could  form  the  basis  for  publishabk* 
research  on  internal  waves.  In  particular,  the  nature  of  the  large-o  (or  k  -  limit  couid  be 
clarified  using  ideas  from  Appendix  (>.  This  applies  not  oiil>  to  the  dispersion  relation,  it 
is  possible  to  use  the  method  of  Appendix  (i  to  determine  ail  om  output  {larameter.'  in 
the  short-wave  limit,  but  this  has  not  yet  betui  done.  There  is  also  a  numerical  eiemeni 
to  this:  our  program  has  difficulty  in  distinguishing  modes  in  the  shon-wave  hnut  1 
above).  If  it  could  be  guided  by  the  a-^ymptotic  anahuic  work,  it  might  be  atile  to  do  that 
job  more  systematically, 

6.  Graphical  output 

There  has  not  la-vn  time  to  do  as  good  a  display  job  as  we  would  have  liked  on  tin- 
output  data.  Several  options  were  investigated  bnedv.  e.g.  the  Sun  V'lsuin  package.  .Mat  lab 
etc.  We  shoudd  be  able  to  produce  quite  nice  views  of  the  flow  field,  including  "rnb^x-r 
sheet "  plots,  perspective  views  from  any  direr-rion.  etc.  An  iniporranT  riew  of  the  output 
is  the  "cellular"  or  "looped"  structure  in  the  cross-flow  (.r=constant  i  plane,  dismissed  m 
Appendix  11.  We  have  not  had  time  to  prepare  accurate  plots  of  this  structure.  Another 
form  of  output  not  pursued  is  the  actual  wave  displacement .  e.g.  particle  motions. 

7.  Other  output  and  detection  matters.  Indun  t  etfecr>  of  the  inn  rnai  wave 
field  should  be  given  more  time.  These  include  electromagnet ic  eifecis  etf('ct>  on  capilhuT.- 
wave  size  and  location,  and  other  determinants  of  the  radar  return 
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EOT  post-mortem 

"Post  mortem"  is  not  an  accurate  description  of  the  fate  of  this  project .  I  ho[)el  There 
is  much  more  that  can  and  should  be  done. 

In  the  beginning,  everyone  involved  had  a  lot  to  learn  about  wakes  and  interntii  waves 
If  we  knew  more  to  begin  with,  we  could  have  avoided  some  pitfalls.  We  might  in  fact 
have  chosen  to  do  little  of  what  was  eventually  done,  relying  iirstead  for  example  t)n  the 
Canadian  study  (Appendix  S).  which  adopt.?  a  similar  methodology. 

I  believe  that  would  have  been  a  mistake  It  was  necessary  for  our  own  computer 
program  to  be  developed,  in  order  that  we  could  build  up  onr  owti  knowledge  and  nnuitton 
This  has  now  been  done.  It  is  possible  that  our  program  is  better  than  those  of  {previous 
investigators.  It  is  certciinly  comparable,  and  is  producing  results  of  comparable  cpiality. 
This  in  less  than  2  months  at  a  very  small  total  cost. 

Internal  waves  are  not  very  big.  Some  of  us  would  like  to  have  sc'en  waves  with 
velocities  of  the  order  of  metres  per  second  instead  of  millimetres  per  sectmd.  It  is  not  for 
me  to  judge  whether  these  are  measureable  or  detectable  signals.  Certainly  it  appears  that 
detection  must  be  indirect,  and  that  bzis  also  been  appreciated  by  a  number  c.^  previon,- 
investigators. 

.A,t  the  beginning,  it  was  felt  that  more  of  my  time  might  be  spent  on  nori-interna! 
wave  aspects  of  the  wake.  This  has  not  proved  possible  in  a  direct  way.  However,  indirectly 
we  are  covering  some  parts  of  the  general  problem.  In  a  crude  sense,  the  viscous  wake  can 
be  interpreted  as  a  "tail"  added  to  the  stern  of  the  submarine,  in  which  ctise  one  can  get 
information  about  its  effect  by  studying  a  "non- balanced"  source-sink  pair  in  the  Rankine- 
ovoid  model.  Immediate  properties  of  the  actual  turbulent  wake  are  another  matter. 

From  the  personal  point  of  view,  this  has  been  a  most  challenging  and  interesting 
project.  I  thank  David  Cartwright  for  organising  it.  and  my  colleagues  while  working  on 
it.  Graham  Furnell.  Tony  Legg  and  Michciel  Carroll.  I  have  not  been  ea-sy  to  work  with, 
being  ever  conscious  of  the  limited  time  available  for  me  to  contribute.  Now  that  thi.s  time 
is  up,  I  hope  that  the  work  can  be  carried  on  within  DSTO. 
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Final  results. 

This  secnon  of  the  reptort  is  wTirten  .-.eptarately  from  the  rest,  m  p;in  berausp 
minute  interruptions  to  the  eompuring  system  that  p>n>vented  us  mco’-;  .  .  ip..  mn)  the 
main  report,  and  in  parr  because  the  same  interruptions  prevented  us  from  obtaining  iin 
time)  as  high  a  quality  of  final  output  as  we  would  have  preferred.  It  is  piamied  that  a 
separate  report  will  be  produced  by  Graham  Furnell  when  a  compleo'  .<et  of  fully-ch(*<:k(*d 
final  results  are  obtained. 

By  ''hnal"  output,  we  mean  actual  computations  of  the  flow  field  [u.  c.  a)  conifKiiieiits. 
as  functions  of  the  spatial  co-ordinates  (x. ,(/.  c ).  This  "only"  requires  nurnenral  integration 
of  the  integrals  appearing  in  .\ppendix  12,  since  now  are  coiihdeni  ui  our  atnhry  n:) 
compute  accurately  the  generation  amplitudes  .4  and  B. 

However,  there  are  pitfalls  in  that  numerical  inregration  task.  The  mam  dithculn-  is 
that  when  x  and  y  are  large,  the  integrand  of  this  6/-integratioii  is  rapidly  varying.  This 
is  no  surprise,  of  course:  it  is  the  basis  for  Kelvin  s  stationary  plimse  argument.  There  ai'c 
several  methods  to  account  for  it  numerically.  One  (set-  Tuck.  Collins  and  Wells  1971 1  is 
to  mimic  the  method  of  stationary  phase.  This  wa-;  used  by  Madiira^inghe  in  producing 
Figure  1.2. 

■Another  more  brute-force  method  is  to  use  many  values  of  h  in  the  discretisation. 
However,  we  don't  want  or  need  to  compute  A  and  B  at  man>-  H  (hence  rr  and  k)  values, 
since  these  are  not  very  rapidly  \-arying  functions,  though  e.xpensive  to  compute.  It  is 
only  the  sine  and  cosine  functions  involving  x  and  y  that  are  rapidly  var>  ;:.!i.  Hence  a 
good  procedure,  also  tried  by  Madurasinghe.  is  to  interpolate  within  (sayi  30  computed 
vulues  of  .4.  while  dividing  the  integration  range  of  the  M  integral  into  many  thousands  of 
sub-intervals. 

We  ha%’e  had  time  to  do  neither  of  these  things  yet.  Hence  the  results  attached  are 
only  reliable  for  re.  uvely  small  values  of  x  and  y.  In  particular,  it  is  not  possible  to 
demonstrate  any  asymptotic  decay  as  x  —  or  y  —  tc  unless  this  rapid  oscillation  is 
accounted  for  numerically,  so  our  resirlts  display  spurious  oscillations  that  do  not  decay  at 
infinity. 

The  Figures  attached  show  the  lateral  velocity  vix.yJ))  at  the  free  surface,  as  a 
function  of  y  for  fixed  x.  The>-  were  computed  for  a  submarine  of  length  L  =  100  metres 
and  radius  /?  -  5  metres,  travelling  at  U  =  2.-5  metres /set'ond  at  depth  h  =  35  metres  m 
an  ocean  with  the  stratification  of  Dawson's  report.  These  rt>sults  were  computed  using 
Simpson's  rule  on  either  30  or  60  interwals. 

For  example,  the  first  two  figures  are  at  x  =  500.  which  is  probably  alreadt-  too 
large,  with  30  and  60  f^-points  respectively.  The  30-point  output  is  smoother,  for  reasons 
not  entirely  appreciated  at  this  time,  and  there  is  general  order  of  magnitude  agreement. 
There  is  good  quantitative  agreement  between  these  curves  about  the  major  minimum  of 
about  —0.0007.5  (nearly  a  millimetre  per  .second)  at  about  y  =  25  metres  and  (perhaps)  a 
subsequent  maximum  of  -f-0.0002  at  about  y  =  100  metres,  but  this  quantitative  agreement 
vanishes  for  larger  y.  and  the  resirlts  are  no  longer  to  be  beheved  at  these  y  values. 


Wf  have'  conipured  some  other  oiitimt  at  iiu>re  reasonable  isiiialiem  value'  of  c  rten- 
erating  waves  of  the  order  of  i.tue  millimetre  per  second,  but  at  the  nine  of  writing  have 
not  bmi  able  to  plot  it.  Plots  will  fn'  attached  if  pnalm-ed  before  the  end  of  T!il^  da\ 
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Figure  15.5:  u{x,y,Q)  versus  y  for  a:  =  0  to  100  meters. 
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Velocity  Component,  W  vs  Y  at  X  =  0  ..  100  (Step  =  50)  Z  =  30m 

30  Theta  Values,  Y-Grid  =  Im,  Submarine  Depth  =  35m,  Submarine  Speed  =  2.5  m/s 
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Figure  15.7:  t'(a',2/,0)  versus  y  for  a-  =  0  to  100  meters.  Speed  ('  =  10  ms-\  Note 
smaller  wave!  (and  almost  independent  of  i  -  to  be  expected,  see  formula:  more  and 
more  diverging  waves  as  U  -*  oo). 
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Velocity  Component,  V  vs  Y  at  X  =  0 ..  100  (Step  =  25)  Z  =  Surface 

30  Theta  Values,  Y-Grid  =  5m,  Submarine  Depth  =  35m,  Submarine  Speed  =  1 .5  m/s 
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Figure  15.S:  u(x,y,0)  versus  1/  for  x  =  0  to  100  meters.  Speed  U  =  I..*}  ms-\  Note  biggir 
wave!  (“Critical  speed’’  where  k  =  <ti  ~  4.7  is  U  =  1.5!)  i.e.  this  is  when  =  0. 
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